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ABSTRACT 


This  study  is  concerned  with  optimal  policies  for  a  sequence 
of  discrete  capacity  expansions  in  response  to  expected 
growth  in  demand  over  an  arbitrary  planning  interval. 

Demand  growth  is  an  arbitrary,  differentiable,  strictly 
increasing  function  of  time.  We  assume  economies  of  scale 
for  investment  costs.  Time  intervals  between  capacity 
expansions,  as  well  as  the  magnitudes  of  the  expansions,  are 
continuous  variables.  Discounting  is  continuous. 

The  model  used  for  a  single  capacity  expansion  is  a  member 
of  an  isomorphic  family.  The  first  model  includes  invest¬ 
ment  and  shortage  costs;  the  second  includes  investment, 
shortage  and  production  costs;  and,  the  last  includes 
investment  costs  and  revenue.  The  first  model  was  chosen 
for  analysis. 

An  asymptotic  stationary  model  is  derived  for  a  linear  demand 
function.  This  model  assumes  that  capacity  expansions  are 
of  equal  size. 

Optimal  policies  have  the  following  property:  either  an 
expansion  of  size  zero  is  optimal  or  capacity  must  be 
expanded  to  equal  or  exceed  demand  at  the  time  of  an 
expansion. 

By  detailed  analysis  of  a  single  expansion  problem,  an 
algorithm  is  developed  for  determining  both  the  optimal  size 
and  time  of  installation  for  the  expansion.  For  a  sequence 
of  expansions  over  some  planning  horizon,  a  similar  calcula¬ 
tion  is  intractable.  However,  by  taking  discrete  values  for 
the  size  of  expansion,  which  is  reasonable,  it  is  possible 
to  use  part  of  the  algorithm  for  the  single  expansion  problem 
to  determine  optimal  sizes,  times  and  number  of  expansions 
over  any  planning  horizon. 

Computations  of  optimal  sizes,  times  and  number  of  expansions 
are  made  for  various  demand  curves.  These  computations  show 
that  the  size  and  time  of  the  first  expansion  (as  a  function 
of  the  planning  interval)  converge  to  fixed  values.  In  the 
case  of  linear  demand  laws,  it  is  shown  that  the  general 
treatment  used  converges,  as  the  planning  interval  tends  to 
infinity,  to  the  solution  for  a  stationary  model. 

It  is  concluded  that  the  models  analyzed  are  susceptible  to 
practical  application,  because  the  solution  for  the  first 
expansion  is  not  sensitive  to  the  demand  law  used. 


GLOSSARY  OF  FREQUENTLY  USED  SYMBOLS 


D(y)  **  Demand  function 

F(  )  -  Inverse  function  of  D 

I  »  Length  of  planning  interval 

v  -  Initial  capacity  at  beginning  of  planning  interval 
x  ■  Size  of  a  capacity  expansion 
t  «*  Time  at  which  a  capacity  expansion  occurs 
n  *  Number  of  expansions 

k  *  Investment  cost  constant 

a  ■  Economies  of  scale  constant  for  investment  costs 

p  **  Penalty  cost  constant 

r  =  Interest  factor  for  present  value  computations 
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CHAPTER  I 


INTRODUCTION 

This  paper  is  concerned  with  optimal  policies  for  a  sequence  of  discrete 
capacity  expansions  in  response  to  expected  growth  in  demand  over  an 
arbitrary  time  interval.  Demand  refers  to  a  single  commodity.  For  each 
expansion,  there  are  two  decision  variables:  the  size  of  the  expansion  and 
the  time  at  which  it  occurs.  Capacity  expansion  over  time  is  a  step  function. 
Discounting  is  continuous.  The  problem  is  analogous  to  a  class  of  inventory 
problems  in  which  order  quantities  are  substituted  for  capacity  increases. 

Several  investigators  have  studied  aspects  of  the  general  problem  of 
capacity  increases.  Since  problems  associated  with  continuous  expansion  have 
been  widely  studied,  we  mention  here  only  a  few  references.  Arrow,  Beckman 
and  Karlin  [  1  ]  studied  a  model  in  which  capacity  may  change  continuously 
(positively  or  negatively)  over  time  with  an  upper  bound  on  the  rate  of 
change  at  any  time.  They  assumed  constant  returns  to  scale  for  all  operating 
and  investment  costs.  Bradley  [  2  ]  studied  a  multi-commodity  continuous 
expansion  problem. 

Hinomoto  [  7 ]  studied  some  effects  of  decreasing  investment  and 
operations  costs  over  time  which  result  from  technological  improvements. 

Most  of  the  other  investigations  assume:  capacity  expansion  to  be  a 
step  function;  obsolescence  of  plant  is  omitted  from  consideration;  lead  time 
is  not  considered;  economies  of  scale  exist  for  plant  investment  costs. 
Following  is  a  brief  summary  of  the  more  recent  studies. 

Chenery  (  3)  assumed  no  shortages,  a  constant  rate  of  growth  in  demand 
and  equal  spacings  between  a  finite  number  of  capacity  increases. 

Howard  and  Nemhauser  [  8 ]  consider  a  capacity  expansion  model  analogous 
to  an  inventory  model  with  demands  known  and  capacity  changes  occurring  at 
discrete  points  in  time.  They  assume  a  composite  convex  function  for  the 
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combined  shortage  and  excess  capacity  costs,  hence  economies  of  scale  for 
investment  costs  are  not  explicitly  considered.  The  model  approximates  the 
situation  in  which  continuous  capacity  expansions  are  permitted.  Manne  [  9] 
and  [10]  studied  a  stationary  policy  for  a  linear  demand  function  over  an 
infinite  planning  interval.  He  also  considered  a  stationary  model  in  which 
demand  is  a  discrete  random  variable  with  an  infinite  planning  interval. 

Manne  and  Veinott  [10]  studied  a  model  which  assumed  that  capacity 
increases  occur  at  given  points  in  time,  nondecreasing  demand  with  shortage 
defined  as  the  difference  (positive)  between  demand  and  capacity  only  at  the 
time  of  a  capacity  increase. 

Srinivasan  [10]  studied  a  model  which  assumed  no  shortages  and 
exponential  demand  growth  over  an  infinite  planning  interval. 

Erlenkotter  [10]  studied  a  model  which  assumed  an  Infinite  planning 
interval,  constant  growth  rate  of  demand,  and  shortage  defined  as  the 
difference  (positive)  between  demand  and  capacity  only  at  the  time  of  a 
capacity  increase. 

The  U.  N.  Bureau  of  Economic  Affairs  [11]  published  a  study  of  a  single 
expansion  model,  with  piece-wise  linear  demand  growth  without  discounting, 
in  which  an  optimal  policy  is  one  which  minimizes  the  investment  cost  per  unit 
of  output  over  a  finite  planning  interval. 

In  this  study,  we  allow  demand  to  grow  in  an  almost  arbitrary  way  to 
account  for  the  large  variety  of  growth  patterns  which  appear  to  be  typical. 
For  example,  the  world  demand  for  nitrogenous  fertilizer  during  1955—1960  and 
cement  during  1947-1960  grew  at  decreasing  rates  [12].  Hence  a  power 
function  with  an  exponent  less  than  one  is  a  convenient  way  to  describe 
demand  as  a  function  of  time.  A  contrary  condition  prevails  for  glass 
container  demand  in  Central  America  where  demand  in  1959  was  growing 
geometrically  [11].  It  is  not  uncommon  to  observe  demand  growing  at  an 
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increasing  rate  at  first  but  at  a  decreasing  rate  later.  A  representation 
of  such  a  demand  pattern  takes  the  form  of  an  "ogive,"  a  curve  often  used  to 
describe  biological  growth.  Therefore,  this  study  assumes  that  demand  is  a 
strictly  increasing,  differentiable  function  in  order  to  account  for  a  wide 
variety  of  growth  situations.  Monotonicity  is  required  only  for  the  unique¬ 
ness  properties  of  the  solutions.  At  the  beginning  of  the  planning  interval, 
Initial  capacity  may  be  larger  or  smaller  than  demand  and  initial  demand  is 
nonnegative. 

An  important  assumption  of  most  pre-investment  studies  is  that  invest¬ 
ment  costs  are  subject  to  economies  of  scale,  e.g.,  double  the  capacity  at 
less  than  double  the  cost.  This  is  a  wide  spread  phenomena.  Haldi  and 
Whitcomb  [  6 ]  report  on  221  cases  of  various  kinds  of  industrial  equipment  in 
14  countries,  84%  of  which  were  characterized  by  an  economies  of  scale 
parameter  of  "a"  which  is  less  than  one,  with  a  modal  value  of  0.6.  This 
seems  to  support  the  "two-thirds  rule-of-thumb,"  commonly  used  in  engineering 
studies.  In  United  Nations  [ll],  the  values  of  0.75  and  0.6  were  used  for  a 
bottling  plant  and  a  fertilizer  plant,  respectively.  Chilton  [  4]  reviewed 
36  types  of  chemical,  petroleum  and  primary  metal  plants  and  concluded  that  a 
value  of  0.6  is  safe  fo  use  for  pre-design  economic  studies.  In  this  study, 
we  assume  that  a  facility  of  capacity  "x"  is  constructed  (or  expanded)  at  a 

g 

cost  of  kx  where  k>0,0<a<l.  This  assumption  was  made  by 
Manne  [  9]  and  [10],  Srinivasan  [l0]>  Erlenkotter  [10]  and  the  United  Nations 
Bureau  of  Economic  Affairs  [11]. 

Economies  of  scale  for  operating  costs  (excluding  costs  for  capital) 
also  occurs  for  many  industries  but  appears  to  be  less  widespread  and  often 
less  important  than  economies  of  scale  for  investments.  For  many  industries, 
consumption  of  fuel,  electric  power,  raw  materials  are  almost  proportional  to 
production.  Capital  intensive  industries,  such  as  process  industries,  tend 


V 


to  exhibit  lower  labor  requirements  per  unit  produced  as  plant  size  increases; 
however,  for  these  industries  labor  cost  is  often  small  compared  to  capital 
cost.  Accordingly,  primary  emphasis  in  this  paper  is  given  to  economies  of 
scale  for  investment  costs. 

Since  the  time  intervals  between  expansions  are  continuous  variables,  we 
assume  continuous  compounding  with  an  interest  factor  of  "r"  so,  the 
present  value  of  1  expended  t  years,  hence,  amounts  to  e  rt  . 

The  number  of  increases  in  capacity,  "n"  is  a  parameter  of  significant 
Interest.  For  a  finite  planning  horizon,  as  "n"  becomes  large,  the  size 
of  some,  perhaps  most,  of  the  increases  could  become  small  and  in  any  real 
situation,  it  may  not  be  sensible  to  permit  very  small  Increases.  Accordingly, 
one  could  set  a  lower  bound  on  the  size  of  a  capacity  increase,  however,  we 
omit  this  from  explicit  consideration. 

In  summary,  the  class  of  problems  which  are  considered  herein  are 
characterized  as  follows: 

(i)  demand  per  unit  time  is  a  monotone  Increasing,  continuously 
differentiable  function  of  time. 

(11)  capacity  expansion  is  a  step  function. 

(ill)  shortages  and  excesses  of  capacity  are  permitted. 

(iv)  the  time  intervals  between  capacity  expansions  as  well  as  the 
magnitudes  of  increases  are  continuous  variables. 

(v)  discounting  is  continuous. 

(vi)  the  planning  interval  is  an  arbitrary  finite  length  of  time. 

(An  Infinite  planning  interval  is  also  included  for  a  linear 
demand  function.) 

(vii)  investment  costs  are  subject  to  economies  of  scale. 

(viil)  cost  components  for  optimization  are  the  present  values  of 

Investment  cost  and  shortage  cost  or  production  cost  or  net 


revenue . 
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In  Chapter  II,  we  consider  a  family  of  single  expansion  models.  The 
first  model  is  concerned  with  optimal  policies  for  investment  and  shortage 
costs — shortage  cost  arising  as  a  result  of  penalties  for  supplying  demand 
via  some  expensive  alternative  means  such  as  renting  operating  facilities, 
or  importing  commodities.  The  second  model  does  not  penalize  shortages 
directly  but  instead  is  concerned  with  optimal  policies  for  investment  cost 
and  revenue  produced  by  the  installed  capacity.  The  third  model  considers 
investment  costs,  shortage  costs  and  production  costs.  These  models  are 
shown  to  be  isomorphic  to  one  another;  consequently,  only  one  is  used  for 
subsequent  analysis  and  computation. 

Chapter  III  considers  in  detail  the  single  expansion  model.  The  optimal 
time  of  expansion  is  considered  as  a  function  of  the  size  of  the  expansion. 
Then,  the  optimal  xize  of  expansion  is  considered.  Finally,  an  algorithm 
is  given  for  finding  the  optimal  policy. 

Chapter  IV  considers  a  sequence  of  expansions  over  time.  A  discrete 
dynamic  programming  algorithm  is  developed. 

Chapter  V  considers  an  asymptotic  stationary  model  for  a  linear 
demand  function.  It  is  shown  to  be  a  limiting  case  of  a  sequential 
expansion  model  of  Chapter  IV.  The  resulting  equations  are  the  same  as 
found  by  Manne  [9]  and  Erlenkotter  [10]  using  quite  a  different  approach. 

In  Chapter  VI  we  present  some  computational  results  for  single  and 
sequential  expansion  problems  among  which  is  a  direct  comparison  of  some 
finite  horizon  expansion  policies  with  an  asymptotic  stationary  policy. 
Although  the  number  of  computations  is  limited,  some  interesting 
phenomena  appear.  The  most  interesting  is  that  the  size  of  the  first 
expansion  (as  a  function  of  the  planning  interval)  and  the  time  of 
installing  the  expansion  converge  to  fixed  values.  Moreover,  for  a 
given  set  of  parameters,  the  values  are  the  same  for  a  family  of 
demand  functions. 
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CHAPTER  II 

A  FAMILY  OF  SINGLE  EXPANSION  MODELS 

In  this  chapter  we  consider  a  family  of  models  which  characterize 
common  investment  situations.  We  assume  shortages  and  excesses  of  capacity 
are  permitted.  For  all  models  a  single  expansion  of  size  "x"  occurs  at 
time  "t"  ,  t  e  [0,1]  and  the  present  value  of  the  investment  cost  is 
I(x,t)  where: 

I(x,t)  ■  kxae  rt 

where  "k"  >  0  ,  is  the  cost  per  unit  of  capacity,  "a"  is  a  dimensionless 
scale  factor,  0  <  a  <  1  ,  and  "r"  is  the  dimensionless  interest  factor. 

Let  "v"  be  the  initial  capacity  at  time  zero  and  D(t)  be  the 
demand  at  time  t  .  Then  D(0)  is  the  initial  demand  at  time  zero.  We 
permit 

D(0)  ■  v  . 

Let  F(*)  be  the  inverse  function  of  D(t)  defined  as  follows: 

F(v)  -  minty  |  D(y)  >  v  ,  y  e  [0,1]}  . 

c 

From  this  it  follows  that: 

F(D(y))  -  y 
D(F(v  +  x))  *  v  +  x 

F(v)  -  0  if  v  <  D(0)  . 


Figure  1  shows  these  relationships  for  one  of  several  possible  configurations. 
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Define : 


L 

U(x,t)  «=  j  min  (D(y),v)e  ^dy 
0 

I 

+  J  min  (D(y),v  +  x)e  rydy  . 


The  general  form  for  the  present  value  of  total  cost  for  the  family  of  single 
expansion  models  is: 


J(x,t)  «  b^I (x ,t)  +  b2U(x,t)  +  b^ 


where  »  b2  »  b3  are  constants.  We  will  illustrate  three  models  which 

are  useful  in  characterizing  investment  problems  and  which  differ  from  one 
another  only  in  the  choice  of  parameters. 

MODEL  ONE 

Assume  capacity  shortage  at  any  time  carries  a  penalty  of  $p  per 
unit  short.  In  Figure  1,  the  shortage  is  the  shaded  area  between  the  demand 
and  capacity  curves.  Let  T(x,t)  be  the  present  value  of  the  total  costs 
for  investment  and  shortages.  Then: 


T(x,t) 


I (x ,  t)  -  pU(x,t)  +  p 


I 

J  D(y)e"rydy 
0 


hence 


I 

bx  -  1  ,  b2  -  -p  ,  b3  -  p  J  D(y)e"rydy  . 

0 


This  first  function  has  an  alternate  form  which  is  more  convenient  to  use. 
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namely: 


t 

T(x,t)  «  kxae  rt  +  p  J  max  [D(y)  -  v,0]e”rjrdy 


0 

I 


+  p  J  max  (D(y)  -  v  -  x,0]e  ^dy  . 


MODEL  TWO 

Suppose  revenue  is  earned  in  proportion  to  the  output  produced  and  that 
production  is  equal  to  the  demand  or  capacity,  which  ever  is  smaller.  For 
example,  the  cumulative  output  in  Figure  la  corresponds  to  the  shaded  area. 
Assume  production  costs  are  proportional  to  output,  hence  revenue  means 
net  revenue.  Let  $q  be  the  revenue  per  unit  produced  and  R(x,t)  be  the 
present  value  of  the  revenue  less  investment  cost  for  an  expansion  of  size 
x  at  time  t  over  the  planning  period  of  length  I  .  Then: 

R(x,t)  ■  -l(x,t)  +  qU(x,t)  . 

Model  One  and  Model  Two  are  simply  related  as  follows: 

T(x,t)  c  -R(x,t)  +  constant  , 

whenever  p  *  q  .  Here  we  may  interpret  shortage  penalty  as  a  loss  in 
revenue. 

MODEL  THREE 

Suppose  that  in  addition  to  investment  cost  and  shortage  cost  of  Model 
One,  we  also  incur  a  production  cost  of  Sc  per  unit  produced.  Assume  as 
before  that  the  amount  produced  is  the  minimum  of  capacity  or  demand.  Let 
P(x,t)  be  the  present  va. -e  of  investment,  shortage  and  production  costs. 
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Then: 


P(x,t) 


I(x,t) 


(p  -  c)’«J(x,t)  +  p 


I 

j"  D(y)e"rydy  . 
0 


This  is  of  the  same  general  form  as  J(x,t)  where 

I 

bi "  1  *  b2  “  "(p  "  c)  •  b3 “ p  j  DMe~Ty*y  • 

0 


SUMMARY : 

We  have  shown  that  three  different  models  of  investment  situations 
belong  to  the  general  family  of  models  defined  by  J(x,t)  and  that  they 
differ  from  each  other  only  in  the  choice  of  parameter  values.  The  properties 
of  this  family  may  be  derived  from  or  discussed  in  terms  of  any  model  in 
the  family,  in  particular,  the  first  model  with  cost  function  T(x,t)  . 
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CHAPTER  III 

ANALYSIS  OF  THE  SINGLE  CAPACITY  EXPANSION  PROBLEM 

COST  AS  A  FUNCTION  OF  "t" 

Consider  the  cost  function  T(x,t)  .  We  will  use  the  following 
abbreviated  notation.  Let  T(t  |  x)  be  the  cost  function  for  a  given  value 
of  x  and  let 

T'(t  |  x)  T(t,x) 

with  similar  notations  for  other  partial  derivatives.  Recall: 

v,0)e  rydy 

v  -  x,0)e  r^dy  . 

In  general  the  first  partial  derivative  with  respect  to  t  is: 


T'(t  |  x)  -  -rkxae  rt  +  pe  rtmax(D(t)  -  v,0) 

(la) 

-  pe  r  max(D(t)  -  v  -  x,0)  . 

c 

Assume  that  v  <  D(I)  and  that  x  is  restricted  so  that  F(v  +  x)  <  I  . 
Hence  F(v)  <  I  .  The  relative  magnitude  of  the  initial  conditions  v  and 
D(0)  provide  two  cases:  D(0)  <  v  and  D(0)  >  v  .  For  each  case  there 
are  three  intervals  for  t  for  a  given  value  of  x  . 

Case  A:  v  >  D(0) 

(i)  t  c  [ 0 ,  F(v)  ] 

or  D(t)  -  v  <  0 

T'(t  |  x)  -  -rkxae"rt  <  0  . 


Hence : 
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(il)  t  e  (F(v) ,F(v  +  x)  ] 
or  D(t)  -  v  >  0 

and  D(t)  -  v  -  x  <  0 

Hence : 

(iii)  t  c  (F(v  +  x)  ,l] 

or  D(t)  -  v  -  x  >  0 

Hence : 


T’(t  |  x)  -  e_rt[-rkxa  +  p(D(t)  -.v)J 


T'(t  (  x)  -  e~rt[-rkxa  +  px]  . 


Remark : 


T'(t  |  x)  is.  continuous  for  all  t  e  [0,1]  because 


T'  (t  |  x) 
T'(t  |  x) 


-  -rk*VrF(v) 


t-F(v) 


t«F(v+x) 


(-rkx‘  +  px)e_rF<vrt)  . 


Consider  two  intervals  for  x  : 


Al:  0  <  x  < 

sr  m 


(f) 


1_ 

1-a 


A2: 


(f) 


1_ 

1-a 


<  x  <  D (I )  -  v  . 


1_ 

1-a 


Define 


•  -  m  ■ 


Case  Al: 


<  X  <  W  ■>  DX  < 


px'‘  <  rk  , 


or 


px  <  rkx  . 


Hence 
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T'(t  |  x)  <  0  for  t  c  [F(v  +  x) ,1)  . 

Now  [-rkxa  +  p(D(t)  -  v]  is  strictly  increasing  in  t  for  all 
t  e  [F(v),F(v  +  x)]  hence 

T'(t  |  x)  <0  for  all  t  e  [F(v) ,F(v  +  x)]  . 

* 

Let  t  (x)  be  the  best  choice  of  t  for  a  given  x  .  Then: 


t  (x)  ■  I 


for  0  <  x  < 


1_ 

1-a 


w  . 


Case  A2: 


The  inequality  of  Case  A2  implies: 


px  >  rkx  . 


Hence 


T'(t  |  x)  >0  for  t  e  [F(v  +  x),I]  . 


But  from  (i)  we  have 


T'(t  lx)  <0  for  t  e  [0,F(v)]  . 


Hence  T'(t  |  x)  changes  sign  from  negative  to  positive  for 
t  c  (F(v),F(v  +  x))  and  since  [-rkx  +  p(D(t)  -  v)]  is  strictly  increasing 
in  t  it  follows  that  there  exists  a  unique  t  (x)  satisfying 


D(F(v))<  D(t*(x)  )■  ^yjxa  +  v  <  D(F(v  +  x)) 


and 


'<>)  ■  f(^  ♦  v) 


for  w  <  x  <  D(I)  -  v  . 
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Summary  of  Case  A:  v  >  D(0) 


t*(x) 


'I  for  0  <  x  <  min  {w,D(I)  -  v} 


Xa  +  vj  for  w  <  x  <  D(I)  -  v 


Note:  for  t  (x)  <  I  we  have: 


(2)  F(v  +  w)  <  F(v  +  w*  axa)  <  F(v  +  x)  . 


Case  B:  0  <  v  <  D(0) 


for 


1 

Denote  z  ■  (D(0)  -  v)J  .  Here  we  must  consider  several  intervals 

x  depending  on  the  relative  size  of  w  and  D(0)  -  v  as  follows. 


B(l)a:  0  <  x  <  w  <  D(0)  -  v 

B(l)b:  0  <  w  <  x  <  D(0)  -  v 

B(l)c:  0  <  w  <  D(0)  -  v  <  x  <  z 

B(l)d:  0  <  w  <  D(0)  -  v  <  z  <  x  <  D(I)  -  v 

B(2)a:  0  <  x  <  D(0)  -  v  <  w 

B(2)b:  D(0)  -  v  <  x  <  w 

B(2)c:  D(0)  -  v  <  w  <  x  <  D(I)  -  v 

Note:  Cases  B(l)a  through  B(l)d  assume 

w  <  D(0)  -  v 


whereas  Cases  B(2)a  through  B(2)c  assume 

w  >  D(0)  -  v  . 

Note:  0  <  v  <  D(0)  ->  F(v)  -  0  hence  for  B(l)a  and‘B(l)b: 


X 
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* 


x  <  D(0)  -  v  ■> 


D(t)  -  v  -  x  >  0 
(D(t)  -  v  >  0 


V  t  e  [0,1] 


hence  T'(t  j  x)  -  e”rt[-rkxa  +  px]  for  all  t  in  the  interval  t  e  [0,1] 
This  single  form  for  t'(t  |  x)  simplifies  the  analysis  for  Cases  B(l)a 
and  B(l)b. 


Case  B(l)a; 


x  <  w  ■£>  px  <  rk 


<0  px  <  rkx 


■>  -rkxa  +  px  <  0 


thus  Equation  (la)  is  nonpositive,  i.e., 


T’(t  j  x)  <  0  V  t  e  [0,1] 


hence 


t  (x)  ■  I 


for  0  <  x  <  w  . 


Case  B(l)b: 


w  <  x  <  D(0)  -  v  . 


q 

This  implies  that  -rk  +  px  >  0  or 


-rkx8  +  px  >  0  , 

thus  Equation  (la)  is  positive,  i.e., 

T’(t  j  x)  >  0  V  t  e  [0,1]  . 

Hence 

t*(x)  °  0 


for  w  <  x  <  D(0)  -  v  . 


Summary  for  Case  B(l)a  and  B(l)b 


I  for  0  <  x  <  w  <  D(0)  -  v 

m  m  m 

0  for  0  <  w  <  x  <  D(0)  -  v 

Case  B(l)c: 

(2a)  0  <  w  <  D(0)  -  v  <  x  <  z  . 

Consider  the  two  intervals  on  t  : 

(i)  t  e  (F(v  +  x),I] 

T’(t  |  x)  -  e'rt(-rkxa  +  px] 

(ii)  t  e  [0,F(v  +  x) ] 

T'(t  j  x)  =  e  rt[-rkxa  +  p(D(t)  -  v) ]  . 
Now  (2a)  implies  -rk  +  px  >  0  or 

-rkx  +  px  >0. 


Hence 

(3)  T'(t  |  x)  >  0  V  t  e  (F(v  +  x),I]  . 

Note,  however,  that  [-rkx  +  p(D(t)  -  v)]  is  strictly  increasing  in  t 
and  -rkx  +  p(D(0)  -  v)  may  be  nonnegative  (when  t  *  0)  .  But  Case  B(l)i 
implies : 


% 


1 

<  D(0)  "  v  =>  ^  < 

=>  (D(0)  -  v)a  <  (D(0)  "  v) 

1 

«>  D(0)  -  v  <  ^  (D(0)  ~  v)]  • 

1 

Hence  for  D(0)  -  V  <  X  <  (D(0)  -  v)]‘  ,  Case  B(l>c  means  that 

rkxa  <  p (D(0)  -  v)  •»  -rkxa  +  p(D(0)  -  v)  >  0 

«=>  T'  (t  |  x)  >  0  at  t  -  0 

Since  T' (t  |  x)  is  strictly  increasing  this  implies 

(4)  t* (t  |  x)  >  0  V  t  e  [0 ,F(v  +  x)]  . 

Inequalities  (3)  and  (4)  together  imply  that 

T' (t  |  x)  >  0  V  t  e  10,1]  . 

Hence:  t*(x)  B  0  for  Case  B(l)c. 

Remark : 

Whenever  (D(0)  -  v)]“  i  »<«  “  »  Cssc  B(1)c  Ca"  bC  Writte” 

w  <  D(0)  -  v  <  x  <  D(I)  -  V 
because  we  have  restricted  x  e  [0,D(I)  -  v]  . 

Case  B(l)d: 


0  <  w  <  D(0)  -  v  <  z  <  x  <  D(I)  -  v  . 
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For  the  two  intervals  on  t  : 

(i)  t  e  (F(v  +  x),I] 

T'(t  |  x)  ■  e-rt[-rkxa  +  px] 

(ii)  t  e  [0,F(v  +  x)  ] 

T’(t  j  x)  »  e"rt[-rkxa  +  p(D(t)  -  v)]  . 
Case  B(l)d  implies 


rk  <  px 


1-a 


rkx“  <  px 


-rkxa  +  px  >  0  ■=>  T'(t  |  x)  >  0 


for 


But  Case  B(l)d  also  implies 


A  <D(°> 


v)  <  xa  ->  p(D(0)  -  v)  <  rkxa 


■=>  -rkx  +  p(D(0)  -  v)  <  0 


->  T’(t  I  x) 


<  0  . 


t-0 


Since  [-rkx  +  p(D(t)  -  v)]  is  strictly  increasing  there 

jit  4 

t  (x)  e  (0,F(v  +  x))  such  that  T’(t  (x)  |  x)  *  0  and  t 


D(t*(x) )  -  —  xa  +  v  <  (v  +  x) 
P 


or 


(5) 


* 

t  (x) 


x  +  v)  <  F(v  +  x) 


t  E  [F(v  +  x)  ,1]  . 


exists  a  value 
(x)  satisfies 


for  Case  B(l)d.  Summary  for  w  <  D(0)  -  v  and  Q  <  v  <  D(0)  : 
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I  for  0  <  x  <  w 


t  (x)  =  \0  for  w  <  x  <  min  {z,D(I)  -  v} 

F(v  +  V)  for  z  <  x  <  (D(I)  -  v)] 


Case  B(2) a: 


0  <  x  <  D(0)  -  v  <  w 


This  implies  that  -rkx3  +  px  <  0  ,  and  since  T'(t  |  x)  °  e  rt(-rkx 
for  all  t  c  [0,1]  , 

T'(t  |  x)  <0. 

Hence  t  (x)  «=  I  for  Case  B(2)a. 


Case  B(2)b : 


D(0)  -  v  <  x  <  w  . 


Consider  the  intervals  on  t 


(i)  t  c  [0,F(v  +  x)] 


T'(t  |  x)  -  e  rt(-rkx3  +  p(D(t)  -  v)) 


(ii)  t  e  (F(v  +  x) ,1] 


T'(t  |  x)  ■  e  rt(-rkxa  +  px)  . 


Case  B(2)b  implies 


1-a  ,  .  .a 

px  <  rk  ■>  px  <  rkx 


->  -rkx  +  px  <  0 


+  px) 


( 


* 


hence  T'(t  |  x)  <  0  for  t  e  (F(v  +  x),I]  .  But  Case  B(2)b  also  implies 


p(D(0)  -  v)  <  px  *=>  -rkxa  +  p(D(0)  -  v)  <  -rkxa  +  px  <  0 

hence 

T'(t  |  x)  <  0  . 
t*0 

Since  [-rkx  +  p(D(t)  -  v)]  is  increasing  in  t  ,  and  since  T' (t  |  x) 
at  t  -  0  as  well  as  t  e  (F(v  +  x),I)  ,  it  follows  that 

T'(t  |  x)  <  0  V  t  c  [0,1]  . 

A 

Hence  t  (x)  c  1  for  Case  B(2)b. 

Case  B(2)c: 

D(0)  -  v  <  w  <  x  <  D(I)  -  v  , 

or 

£ 

-rkx  +  px  >  0  . 


Hence 


T'(t  |  x)  >  0  V  t  e  (F(v  +  x)  ,1]  . 
Now,  since  x  >  (D(0)  -  v)  , 


T'(t  |  x)  <  (-rk(D(0)  -  v)a  +  p(D(0)  -  v)]e“rt  . 

t=0 

But  Case  B(2)c  implies 

D(0)  -  v  <  w  , 

or 


p (D(0)  -  v)  <  rk(D(0)  -  v)a  . 
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Hence 


T’(t  |  x)  <  0  . 

t*=0 


Since  [-rkxa  +  p(D(t)  -  v)  ]  is  strictly  increasing  in  t  there  exists  a 
value  t*(x)  e  (0,F(v  +  x))  such  that  D(t*(x))  «=  ^^jxa  +  v  <  (v  +  x)  or 


t*(x)  »  F^y  xa  +  vj  <  F(v  +  x) 


for  Case  B(2)c.  Summary  for  B(2)a,  B(2)b,  and  B(2)c,  i.e., 


D(0)  v  <  w  <  D(I)  -  v  and  0  <  v  <  D(0)  : 


* 

t  (x) 


I  for  0  <  x  <  w 


F(y  **  +  v)  for  w  <  X  <  D(l)  -  V 


Summary  for  all  cases: 


(i)  t  (x)  =  I  for: 


0  <  x  <  min  {w,(D(I)  -  v) } 


(ii)  t  (x)  =  0  for: 


rmax  [D(0)  -  v,z] 


min  {w,D(I)  -  v}  c  x  <  min  <max  l  >  ,  D(I)  -  v>  . 

I  (min  [w,(D(I)  -  v) ] )  1 


(iii)  t*(x)  =  f(—  xa  +  v)  for: 


min  \max 


L( 


max  [D(0)  -  v,z) 


min  [w,D(I)  -  v] 


,  D(I)  -  v)  <  x  <  (D (I )  -  v)  . 


These  hold  for  both  initial  conditions,  0  <  v  <  D(0)  and  0  <  D(0) 
and  they  may  be  used  to  set  up  computer  codes. 


Graphs  of  t  (x)  : 

There  are  four  possible  situations  for  t  (x)  as  follows.  See 
Figures  2,  3,  4,  and  5.  (These  are  drawn  left  continuous.) 

(i)  w  >  D(I)  -  v 


t*(x) 


D(I)  -  v 


-*•  x 


FIGURE  2 

(ii)  (D (0)  -  v)  <  w  <  (D(I)  -  v) 


t*(x) 


FIGURE  3 


x-a  a  ,  . 

F(w  x  +  v) 


x=v 


F(v  +  w)  . 


Note : 


F^~  (D(I)  -  v)a  +  vj  <  I  . 

The  strict  inequality  holds  because  otherwise  (D(I)  -  v)  **  w  which 
violates  this  case. 


(iii)  w  <  D(0)  -  v  and  z  <  (D(I)  -  v) 


FIGURE  4 


(iv)  w  <  D(0)  -  v  and  z  >  D(I)  -  v 


t*(x) 


0  w  D(0)  -  v  D(I)  -  v 

FIGURE  5 


Remarks ; 

Recall  from  (2),  (5)  and  (6)  that  whenever  t  (x)  -  F(w*  \a  +  v) 
we  have  that 

t  (x)  <  F(v  +  x) 

regardless  of  initial  conditions  v  >  D(0)  or  v  <  D(0)  .  Hence  the 
graphs  of  these  configurations  are  as  follows. 


Demands , 
Capacity 


FIGURE  6a 
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Demand 
Capacity , 


OPTIMAL  x 


Let 

T (x)  =  T(t  (x) 

,x)  .  There  are 

several 

I: 

w  >  D (I)  -  v 

II: 

D(0)  -  v  <  w 

<  D(I)  -  v 

III 

:  w  <  D(0)  -  v 

and  z  <  D(I)  - 

V 

IV: 

w  <  D(0)  -  v 

and  z  >  D(I)  - 

V 

I: 

w  >  D (I )  -  v  ; 

t  (x)  “  I  . 

T(x)  =  kxae  rI  +  p  j  max  (D(y)  -  v,0)e  rydy 

0 


T'(x)  -  akxa_1e_rI  >0  V  x  e  (0,D(I)  -  v]  . 


Hence  min  T(x)  =  T(0)  . 

Case  II:  D(0)  -  v  <  w  <  D(I)  -  v  . 

* 

Case  11(a):  0  <  x  <  w  ;  t  (x)  ■  I 

- .  — -  — — — —  • 

I 

T(x)  *  kxae  rI  +  p  J"  max  (D(y)  -  v,0)e  rydy 

0 
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T’(x)  -  akxa  rI  >  0  V  x  e  (0,w]  . 


Hence  T(x)  is  increasing  for  all  x  e  (0,w)  . 


Case  11(b):  w  <  x  <  (D(I)  -  v)  ;  t  (x)  ■  F(v  +  w*  axa)  <  F(v  +  x)  . 


t*(x) 


K  f 

T(x)  *»  kxae  rt  ^  +  p  J  max  (D(y)  -  v,0)e  rydy 


p  J"  max  (D(y)  -  v  -  x,0)e  ^dy  . 


t  (x) 


Referring  to  Figure  6,  we  note  that  t  (x)  <  F(v  +  x)  ,  y  e  [0,F(v))  => 
D(y)  -  v  <  0  and  y  e  [0,F(v  +  x)]  ■>  D(y)  -  v  -  x  <  0  .  Hence  for  the 
range  w  <  x  <  D(I)  -  v  we  have: 


* 

t  (x) 


* 

t  (x) 


j  max  (D(y)  -  v,0)e  rydy  =  j  (D(y)  -  v)e  rydy 


j  max  (D(y)  -  v  -  x,0)e  rydy  ■  j  (D(y)  -  v  -  x)e  rydy 

t  (x)  F(v+x) 


Hence  the  cost  function  may  be  written: 


T(x)  -  kxVrt  (X) 


t*(x) 


+  P  j  (D(y)  -  v)e-rydy 
F(v) 


+  P  J  (D(y)  -  v  -  x)e  rydy  . 
F(v+x) 


Hence:  T'(x)  «  akxa_1e”rt  <x)  -  *  (e-rF(vfx>  -  e'rI)  . 
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Denote :  T(w  +  0)  -  TOO  I  !  T<0)  -  TCx)|  .  Then  compare  T(w  +  0) 


with  T(0) ,  where: 


T(0)  —  P  j*  (D(y)  -  v)e  rydy  , 

F(v) 

since  y  c  [0,F(v)]  ->  D(y)  <  v  And 

I 

T(»  +  0)  ■  kvVlF('m'>  +  P  |  <My)  -  v)e'rydy 


e'^dy  , 


F(v+w+0) 


because  at  x 


w  +  0  ,  t*(x  +  0)  n  F(v  +  w  +  0)  .  Hence: 


T(w  +  0)  -  T(0)  -  kw  e 


ae-rF(w+v)  +  £W  _  e-rF(v+w)} 


^kwa  _  Eije-rF(v4-w)  +  Et  ftl  , 


kwa  -  f  , 


hence 


:(w  +  0)  -  T (0)  -  ^  e'rl  >  0  , 


or 
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Case  lib:  Conditions  for  T’(x)  >  0 

to 

Recall  that  for  x  e  (w,D(I)  -  v]  we  found  that 
F(v  +  w)  <  t  (x)  *  F(v  +  axa)  <  F(v  +  x) 

and 

F(v  +  w*  axa)  <  F(v  +  w*  a(D(I)  -  v)a)  <  I  . 

Then  we  can  write 

-rF(v+w)  -rt*(x)  -rF(v+w1_a(D(I)-v) a)  -rl 

6  ?  6  ^  C  ^  c  * 

as 

a-1 

And  since  x  is  decreasing  in  x  , 

£  =  kwa-1  >  kxa_1  >  k(D(I)  -  v)®'1  . 

We  now  use  these  results  to  find  a  lower  bound  on  the  derivative  for 
x  e  (w,D(I)  -  v]  as  follows: 

T'OO  >  akCD(I)  -  v)^le-rF(vfv1-a(D<I)-v)a)  +  £  (e-rl  .  e-tF(v«))  _ 


which  is  nonnegative  if 


*(d(I)  -  v) 


1_ae-rF(v+w1_a(D(I)-v)a)  +  ^-rl  _  e-rF(v+w)^  >  Q 


or 


(7)  I 


1  .  |  -rF(v+w)  /  w  \ 

-  -  t  tn  Le  "  a(wi r=^) 


1-3  -rF(v+w1-a(D(I)-v)a) 
e 


■]. 


in  which  case  T(x)  is  nondecreasing  throughout  the  interval  (w,D(l)  -  v]  . 
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» 


An  example; 

The  importance  of  this  result  may  be  seen  by  means  of  an 
Suppose  D(y)  is  linear  and  units  of  time  are  chosen  so  that 

Then 


F(v  +  w)  =  v  +  w  . 


Suppose  v  =  0  ,  D(0)  =  0,a  =  4,r  =  0.1,k=20,p  =  l 
1 

w  a  3  =  4  .  Then  (7),  the  condition  on  I  becomes: 


I  < 

as 


10  in 


A 


-.2/1 


or 


I  <7.7  units  of  time  approximately  . 

s 

Note:  r  corresponds  to  a  bank  interest  rate  of  6%  ,  i.e., 


1.06  *= 


I 


erydy 


and 


r  ■  .1  approximately  . 


* 


Hence,  for  this  example,  I  <  7.7  units  of  time  mesns  x  ■ 
policy. 

T(x)  is  not  strictly  decreasing  for  all  x  e  (w,D(I)  -  v] 

T  .  .  a-1  -rF(v4v1-a(D(l)-v)a)l 

X’(x)  ■  |_ak(D(I)  -  v)  e  J 

x*D(I) -v 


illustration. 
D(y)  *  y  • 


,  then 


0  is  the  optimal 

because 

>  0  . 


Furthermore,  a  sufficient  condition  for  the  existence  of  relative 


* 
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T(x0) 


kx  e 
o 


*  t*(xo) 

a/"'  6‘.>  +  p  J  (D(y)  -  v)e  rydy 


F(v) 

F(v+x  ) 


(D(y)  -  v  -  x  )e  rydy 


-rt  (x  )  * 

kxae  +  T(t  (0),0)  -  px 

o 


L 

J 


e‘rydy 


F(v+xq) 


F(v+x  ) 


-  P 


I 


(D  (y)  -  v)e_rydy 


t*(xo) 


The  integrand  of  the  last  term  has  known  bounds,  because  for 

t*(*0)  -  F(v  +  w1'***)  i  1  i  F(V  +  *o>  • 

V  +  W1_V  <  D(y)  <  v  +  xo 


and 


w1_axa  <  D(y)  -  v  <  x 


Define  A  -  T(xq)  -  T(0)  .  Then 


I-:-  *  -:fr) '  (•- 


T“a/  -rF(v+x 

rl  0 


-  e 


>) 


/x  \1_V  -rF(vfx  )  )\\ 

u  -e  > 


<  a  , 


-rt  (x  ) 


and 


Then 


and 


|kx  e 
o 


_  -rF(v+x) 
-rl  ^  o 

e  -  e 


+  kxa(e 


-rF(v+x  )  -rt  (x  )] 


-  e 


Hence 


1-a 


-rl 


-rF(v+x 


'•Cfj 


1-a 


-r. 


r--  •  [&r  ■  >] 


-rF(v+x  ) 
e  -> 


(%)“■••" > W  -  J 


-rt  (x  ) 
e  »>  A 


* 

ix  *  x. 


But  if 
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then  the  role  of  xq  Is  uncertain  and  T(XQ)  must  be  compared  directly 
with  T(0)  . 

Case  III:  w  <  D(0)  -  v  and  z  <  D(I)  -  v  . 

(See  Figure.  4).  Consider  four  intervals:  (a)  x  e  [0,w]  , 

(b)  x  e  (w,D(0)  -  v]  ,  (c)  x  e  (D(0)  -  v,z)  ,  (d)  x  e  (z,D(I)  -  v]  . 

* 

Case  111(a):  x  c  [0,w]  ;  t  (x)  c  I 

I 

T(x)  =  kxae  rI  +  p  J  max  (D(y)  -  v,0)e  rydy 

0 


T'(x)  =  akx3  rI  >  0  V  x  e  (0,w] 

T''(x)  =  (a  -  l)akxa  2e  rI  <  0  V  x  e  (0,w]  . 

Here  T(x)  is  concave  increasing  in  this  interval. 

Case  111(b):  x  e  (w,D(0)  -  v]  ;  t  (x)  ■»  0 

Note  that  max  (D(y)  -  v  -  x,0)  *  D(y)  -  v  -  x  ,  and 


I 

T(x)  *  kxa  +  p  J"  [D(y)  -  v  -  x]e  rydy  . 
0 

Hence 


T'(x)  -  akxa_1  +  J  (e'rI  -  1) 


and 
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T''(x)  -  a(a  -  l)kx  <  0  . 


Hence  T(x)  is  concave  for  all  x  e  (w,D(0)  -  v]  ,  but  the  sign  of  T'  (x) 
is  uncertain. 

Note  that 


T'  (x) 


,  a-1  -rl 
akw  e 


T'(x)  **  akw 

x«=w+0 


+  £  (e-rI  -  1)  . 


Hence : 


T'(x)  -  T'(x)  -  akwa-1t 1  -  e“rI]  -  £  (1  -  e 

x=w+0  x=w 


-  ^akwa  1  -  ^j(l  -  e  rl) 


(a  -  l)kwa‘1(l  -  e_rI)  <  0  . 


Therefore  T(x)  is  concave  for  all  x  c  [0,D(0)  -  v]  and  if  relative 
minima  exist  in  the  interval  x  e  (w,z]  they  can  only  occur  in  the  interval 
x  e  (D(0)  -  v,z]  . 


Case  I' 1(c):  x  e  (D(0)  -  vtzj  ;  t  (x)  »  0 


T(x)  ■  kxa  +  p  J  (D(y)  -  v  -  x)e  ^dy 

F(v+x) 

(Note ;  T(x)  is  continuous  at  x  ■  D(0)  -  v  .) 


T'(x)  -  ak*8'1  +£  (e'rI  -  e'rF(vfx); 
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* 


T"(x)  =  a  (a  -  l)kxa“2  +  p(e"rF(v  +  x))(^-  F(v  +  x))  , 

and  the  signs  of  both  T' (x)  and  T''(x)  are  uncertain.  (Note  that 
T' (x)  is  continuous  at  x  =  D(0)  -  v  .) 

Case  111(b):  Conditions  for  T'(x)  >  0  and  T(x)  <  0  for  x  e  (w,z] 

Recall  that  for  Case  Ill,  we  have  that  w  +  v  <  D(0)  and  F(v  +  w  +  0)  e  0  . 
And,  since: 

T-  (x)  =  xkx3’1  +  *  e‘rI  -  £■  e-rF(v+x>  , 
r  r 


we  can  write : 


,  -rF(v+w+0)  -rF(v+x)  ^  -rF(vfz) 

1  *  e  >  e  >  e 


Also  recall  that  for  the  interval  x  e  (w,z] 


p  ,  a-1  ,  a-1  _  ,  a-1 

^  =  kw  >  kx  >  kz 
r  ■ 


Thus  we  may  write  bounds  on  the  derivative  as  follows. 

mi /  \  i  P  /  -rF(v+z)  -rl. 

T'  (x)  <  akw  -  (e  -  e  )  , 

( 

the  right  side  of  this  inequality  is  negative  whenever 


-rl  „  ,  -rF'(v+z)  . 

e  <  (e  -  a)  , 


or 

(8) 


I  > 


(v+z) 


Further,  we  can  write  an  upper  bound  on  T'(x)  as  follows: 
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T'(x)  >  akza_1 


£  ^g-rF(v+w+0) 


or 


T'(x)  >  akza_1  -  ^  (1  -  e"rI) 


because  F(v  +  w  +  0)  =0  .  The  right  side  of  this  inequality  is  positive 
whenever 


> 


or 


(8a) 


If  neither  (8)  nor  (8a)  holds  then  (9)  holds,  where 


and  the  sign  of  T' (x)  for  x  e  (w,z]  is  uncertain,  in  which  case  there 
may  exist  relative  minima  in  the  interval  (D (0)  -  v,z]  .  Note  that  sufficient 
conditions  for  relative  minima  for  x  e  (D(0)  -  v,z]  are  that 


T'(x) 


x-D(0)-v 


<  0  and  T' (x) 


>  0 


x»z 


Suppose  the  first  condition  holds.  Then 


a 


» 


or 


I 


I  <  -  In  I - - - - — \  . 

\e-rF(v fz)  -  a(f  )  / 


Therefore  sufficient  conditions  for  relative  minima  for  x  c  (D(0)  -  v,z] 
are  that 


(9a)  —  £n 


V  '  3(d(0)  -  v) 


r^fj  <1<r"n  VrH"z)  -  <~t~j ' 


(Note  that  whenever  (9a)  holds,  this  implies  that  (9)  holds.) 


Remark;  T' (x)  is  continuous  at  x  =  z  . 


I'  (x)  -  akz3'1  +  £  (e-rF(vtz)  -  B-rI) 

x*=z-0  r 


T'  (x) 


ak^’V"  (x) 


+  *  (e-rF<vfx>  .  e-rI)  . 


x=z+0 


x-z+O 


t*(x) 


F(v  +  w^  aza) 


x«=z+0 


*  F(v  +  w^  aw3  ^"(D(O)  -  v))  =  F(D(0))  =  0  , 


hence 


T'(x)  -  T'(x) 

x-z-0  x-z+O 


Case  III(c):  Compare  T(x)  with  T(0)  for  relative  minima  for 


x  c  (D(0)  -  v,z] 
Let  x  be  the  solution  to 


(e-rF(vkx)  _  e-rl} 
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which  yields  the  smallest  value  of  T(x)  over  the  interval  (w,z] 


T(x)  -  T(x) 


and  T(0)  -  T(x) 


x=x 


x-0 


and 


A  *  T(x)  -  T(0) 


Then 


l 

A  -  kx“  +  p  I  (D(y)  -  v  -  x)e-rydy  -  p  I  (D(y)  -  v)e"rydy 

F(vfx)  0 


J. 

J 


or 


F(v+x) 


A  «  wa  _  £  -rF(vhc)  -rl.  f  -ry 

Alcx  rx(e  -e)-pl  (D(y)  -  v)e  ydy  . 

0 


This  can  be  simplified  by  noting  that  x  satisfies 


ra  -  £  w0“rF(vFx)  _  „~rl, 
r 


akx  -  *  x (e  -  e  “)  . 


Then, 

F(v+x) 

A  -  (1  -  a)kxa  -  P  J  (D(y)  -  v)e_rydy  . 

0 

The  integrand  of  the  last  term  is  bounded  as  follows  i 


0  <  y  <  F(v  +  x)  , 
D(0)  <  D(y)  <  v  +  x  , 
D(0)  -  v  <  D(y)  -  v  <  x 


Denote : 
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Substituting  D(0)  -  v  ,  the  lower  bound,  we  can  write  a  bound  on  A  as 
follows : 


4  <  (1  -  a)kia  -£  <D<0>  -  v)(l  -  *-rF(v+x))  . 


The  right  side  of  this  inequality  is  negative  whenever 


(1  -  a)*1-8,;*  <  (D(0)  -  v)(l  -  e-rF(vtx))  , 


or 


(10a) 


€ 


-rF(v+x), 


(1  -  a)  7  <  (1  -  e  , 


because 


Hence: 


(D(0)  -  v)  *  z^1  3 


(1  -  a) 


(t  r< 


(1  .  ^  T('j  .  T(0)  . 


Substituting  x  ,  the  upper  bound,  we  have: 


4  >  fai8  -  *  - 


F(v+x) 


-rlv 
e  ) 


-  px 


J 


e-rydy  , 


or 


i  "a  p  \  -rF(v+x)  -rl,  £  in  -rF(v+x). 

>  kx  -  x(e  -  e  )~  x(l  -  e  )  , 


or 


:  k*a[F  -  fer<F  -  8'rI)] 


41 


The  right  side  of  this  inequality  is  positive  whenever: 


(§r«  -  '-ri> . 


rl 

e  < 


-0 


a-1  I  ’ 


I  <  —  in 


1  -  0 


a-l  • 


Hence: 


(10b) 


I  <  -  in 


r  -0 


->  T (x)  >  T (0)  . 


If  neither  inequality  (10a)  nor  (10b)  holds  then  a  direct  comparison  of 
T(x)  with  T(0)  must  be  made. 

Case  111(d);  x  e  (ztD(I)  -  v]  ;  t  (x)  *  F(v  +  w*  \a)  <  F(v  +  x)  . 


TOO  -  kxVrt  (X) 


+  p  J*  max  [D(y)  -  v  -  x,0]e  rydy 
t  (x) 


* 

t  (x) 


max  [D(y)  -  v,0]e  rydy  . 


F(vfx) 


T'  Cx)  -  akx8-Vrt‘(x)  -  £  (e-rF<v+,l)  -  e‘rI)  . 

r 

This  is  the  same  expression  as  Case  II  for  the  range 

w  <  x  <  D(I)  -  v  ,  t  (x)  -  F(v  +  w1  V)  . 

Note  that  T'(x)  is  continuous  at  x  ■  z  because  F(v  +  w^ 
Case  111(d):  Conditions  for  T’(x)  >  0  for  x  e  (z,D(I)  - 

Note : 

-rF(v+z)  -rt*(x)  -rFCv+w1  a(D(I)-v)a) 
e  >  e  >  e  ■ 

Also: 

kza  ^  >  kxa  ^  >  k(D(I)  -  v)8  ^  . 

» 

Therefore  we  may  write  a  bound  on  the  derivative  as  follows. 
T'(x)  >  akUHI)  -  v)»-VrP('Mvl"a(D(I,'v)a)  (e‘rI  - 


which  is  nonnegative  if 


(ID 


(\l-a 
D(I)  -  v) 


(e-rF(^)  -  e-rI)  _ 

-rFCv+w1  a(D(I)-v)a) 
e 


aza)  -  F(D(0) ) 


-rF(v+z)j 
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Hence,  if  (11)  holds,  T(x)  is  nondecreasing  throughout  the  interval 
(z,D(I)  -  v]  .  This  is  the  same  as  we  found  for  Case  II,  (see  (7)),  except 
that  F(v  +  z)  has  replaced  F(v  +  w)  .  When  this  inequality  is  satisfied, 
relative  minima  do  not  exist  in  the  interval  x  e  (z,D(I)  -  v]  ,  although 
there  may  exist  relative  minima  in  the  interval  x  c  (w,z]  .  Further,  T(x) 
cannot  be  strictly  decreasing  in  the  interval  (z,D(I)  -  v]  because 


T'(x)  ak(DCI)  -  ^-V^W^CDCn-v)8)  >  o  . 

x=D(I)-v 

If  (11)  does  not  hold,  there  may  be  relative  minima  for  x  c  (z,D(I)  -  v)  . 

A  sufficient  condition  for  the  existence  of  relative  minima  for 
x  f.  (z  ,D(1)  -  v)  is 


T'(x) 


=  akz 


a-1 


x=z 


-  ^  (e 
r 


-rF(v+z) 


-rl. 
-  e  ) 


<  0 


or 


<  (e-^O 


-rl 

e 


) 


9 


or 


(lla) 


(Note  that  whenever  (8)  holds,  T'(x)  <  0  for  all  x  e  (w,z]  ,  which 

implies  that  (lla)  holds.)  Let  x  be  the  solution  to 

o 

ae-rF(v4w1_axa)|wj1_a  =  (e-rF(v+x)  _  Q-rI) 
which  yields  the  smallest  value  of  T(x)  over  the  interval  (z,D(I)  -  v]  . 
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Except  for  the  range  on  "x"  ,  this  is  the  same  result  as  we  found  for 
Case  II  where  x  e  (w,D(I)  -  v]  ;  whereas,  for  this  case  x  e  (z,D(I)  -  v]  . 
Accordingly  we  can  write: 
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‘(d(I)  -  v) 


1-3  -rF(v+w1_a(D(I)-v)a)  ^ 

P  > 


which  is  impossible. 

Case  III  Summary: 

Recall  inequality  (8)  ,  (8a) ,  (9) ,  (9a)  ,  (11)  and  the  implications  when 
each  holds: 


j(8)  =  I  >  i  in  (^VFTi2) 


-  a  i 


((8a) : 


=>  T'(x)  <0  V  x  e  (w,z]  , 


0  V  x  e  (w,z]  , 


=>  the  sign  of 


T*  (x)  is  uncertain  V  x  e  (w,z]  , 


J(9a)  :  ~  In 


V  '  3(d(0)  -  v) 


1  - — \  <  1  <  —  £n 

1-a  1  r 


1 


-rF(v+z)  /w 


/  \l“a 

■•&)  ) 


T’(x) 


<  0  and  T'(x) 


x-D(0)-v 


>  0  , 


X“Z 


((11) :  a 


(ouf^) 


1-a  /  -rF(vfz)  -rl 

e  -  e 

> 


-rF(v+w^  a(D(I)-v)a) 


=>  T'(x)  >  0 


Vx  e  (z,D(I)  -  v]  . 


We  conclude  that: 


(a)  If  both  (11)  and  (8a)  hold,  there  are  no  relative  minima  anywhere 

* 

in  the  interval  (0,D(I)  -  v]  and  hence  x  ■  0  . 
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(b)  If  (11)  and  (9)  hold,  there  may  exist  relative  minima  in 

A 

(D(o)  -v,  z],  and  only  in  this  interval.  And,  for  x 
yielding  the  smallest  value  of  T(x)  in  this  interval: 

x*  -  1  If  U-a)  (-f  )S  <  (l  -e-F(v  +  *>  )  . 

*  ^ 

x  =  0  if  I  <  -  In 

“=  r 

and  T(x)  must  be  compared  directly  with  T(0)  if  neither  of 
the  above  tests  is  satisfied. 

(c)  If  (9)  holds  but  C 11 )  does  not,  there  may  be  relative  minima 
in  either  or  both  intervals: 

(D(O)-v,  z],  (z,  D(I)-v).  And,  if  they  exist  in  both 
intervals,  T(0)  must  be  compared  with  T(x  )  and  T(x). 


* 
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(d)  If  (ID  and  (9a)  hold,  then,  there  exist  relative  minima  in 

A 

(D(0)  -  v,z]  and  only  in  this  interval.  And,  for  x  yielding 
the  smallest  value  of  T(x)  in  this  interval: 


if  (1  -  a) 


(i[<  (i  .  c-lF(vfX>)  , 


* 

X  = 


0 


if 


» 


and  T (x)  must  be  compared  directly  with  T(0)  if  neither  of 
the  above  tests  is  satisfied. 

(e)  Relative  minima  exist  exclusively  in  (z,D(I)  -  v)  only  if 
T'  (x)  <  0  for  some  x  c  (z,D(I)  -  v)  and  either  T  (x)  <0 
for  all  x  e  (D(0)  -  v,z]  or  T'(x)  >  0  for  all  x  e  (D(0)  - 
which  requires  that  (11)  does  not  hold  and  (8)  or  (8a)  holds, 
these  conditions  hold, 


* 

x  = 


if  I  >  F(v  +  xq)  +  -  in 


/  \1_a 

Ht) 


and 


/  ^  1-a  a\  1 

o  1£  I  <  F(v  +  w  +  r 


£n 


[“©  J 


But  if 


F^v  +  w1  V)  +  -^  in 


!-a  ,  <  1  -  F(V  +  Xo) 


v,z]  , 
If 
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the  roJc  of  xq  is  uncertain  and  T(xo)  must  be  compared  directly 
with  T(0)  . 

(f)  Otherwise,  (11)  does  not  hold  but  (9a)  does  hold  and,  relative 
minima  may  exist  in  either  or  both  intervals,  (D(0)  -  v,z]  , 

(z,D(I)  -  v)  and  if  they  exist  in  both,  T(0)  must  be  compared 
with  T(xo)  and  T(x)  ,  by  using  the  tests  above  or  by  direct 
compar ison. 

Case  IV:  w  <  D(0)  -  v  and  z  >  D(I)  -  v 


(See  Figure  5).  This  case  is  similar  to  Case  III  for  the  range 
x  e  (w,z]  .  For  Case  IV  we  found  that 


* 

t  (x) 


I  for  x  e  [0,w] 


/ 

t 


0  for  x  e  (w,D(I)  -  v] 


Accordingly,  most  of  the  required  analysis  has  been  done  and  we  may  summarize 
the  conclusions  from  Case  III  which  apply  to  Case  IV  as  follows: 


a.  T(x)  is  strictly  concave,  increasing  for  all  x  in  the  interval 
[0,w]  . 

b.  T(x)  is  strictly  concave  for  all  x  e  (0,D(0)  -  v]  and  T'(x) 
is  continuous  at  x  =  D(0)  -  v  . 

c.  Since  there  are  no  relative  minima  for  x  e  (w,D(0)  -  v]  ,  we 

need  only  consider  the  range  x  e  (D(0)  -  v  ,  D(I)  -  v)  ,  for  which 

* 

t  (x)  =  0  . 

Case  IV(a) :  x  e  (D(0)  -  v,D(I)  -  v]  ;  t  (x)  =  0 


We  can  write, 
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» 


T(x) 


I 

kx3  +  p  J* 

F(v+x) 


(D(y;  -  v  -  x)e_rydy  , 


and 


T'  (x) 


akxa_1  -  £  (e"rF(v+x) 
r 


-rl 

e 


)  • 


Now  at  the  end  of  the  interval, 


T'(x) 


x=D(I)-v 


ak(D(I)  -  v) a_1  >  o 


as  in  Case  III,  and  T(x)  cannot  be  decreasing  throughout  the  entire 
interval  (D(0)  -  v,D(I)  -  v]  . 

The  existence  of  relative  minima  in  this  interval  requires  that 
T'(x)  <  0  for  some  x  c  (D(0)  -  v,D(I)  -  v]  .  As  before,  we  look  for 
conditions  which  rule  out  such  a  possibility,  namely,  conditions  for  which 
T’ (x)  >  0  for  all  x  t  (D (0)  -  v,D(I)  -  v]  .  As  before  we  note  that  for 
x  e  (D(0)  -  v,D(I)  -  v]  : 


1  >  c-rF(W-x)  £  e-rl 


> 


and 


k(D(0)  -  v)3"1  >  kx3'1  >  k(D(I)  -  v)3_1  . 

Hence  a  lower  bound  on  the  derivative  may  be  written  as  follows: 

T'(x)  >  ak(D(I)  -  v)3-1  -  ^  (1  -  e~rI)  . 
a  r 


The  right  side  of  this  inequality  is  positive  or  zero  whenever 
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a(D(I)  -  v)3"1*1"3  >  (1  -  e‘rI)  , 


or  as  a  condition  on  I  : 


(12) 


!  <  —  in  / — -\ 

r  \  -  w^n 


Hence  we  conclude  that 


I  <  -  in 
s  r 


4 


I  w  \  a(| 

/ 


=>  T'(x)  >  0  Vx  e  (D(0)  -  v,D(I)  -  v] 


Suppose  (12)  holds.  Then,  since 


(D(0)  -  v)3"1  >  (D(I)  -  v)3"1 


we  can  write 


T’(x) 


x=D(0)-v 


ak(D(0)  -  v) 3  1  -  (1  -  e  rI) 


which  is  strictly  positive.  But  T(x)  is  concave  for  all  x  e  (w.D(O)  -  v]  ; 
hence,  whenever  (12)  holds,  T(x)  is  strictly  increasing  for  all 
x  t  (w,D(0)  -  v]  .  Since  T(x)  is  also  concave  increasing  for  all 
x  e  (0,w]  ,  it  follows  that  whenever  (12)  holds,  T(x)  is  nondecreasing 
over  the  entire  interval  (0,D(I)  -  v]  ,  We  conclude  that: 


and  when  this  inequality  does  not  hold,  there  may  exist  relative  minima  in 
the  open  interval  x  e  (D(0)  -  v,D(I)  -  v)  .  A  sufficient  condition  for  the 


existence  of  relative  minima  in  (D(0)  -  v,D(I)  -  v]  is  T' (x) 


<  0 


x=D(0)-v 


* 
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or 


or 


(12a) 


1 


Ut  x  bo  the  solution  to 


which  yield,  the  mlLM  **!»•  of  T(»)  owr  th.  op.i>  lnt.ry.l 

(0(0)  -  v.D(I>  -  v)  .  But  ;  i.  not  th.  o?tl«l  v.lu.  unlc.  T(«)  <  T(0) 

where 

T(i)  -  T(x) I  . 

|x«x 

T(0)  ■  T(x)  I 

|x-0 


Denote: 


A  ■  T(x)  -  T(0) 


A 

*e  p  * /  -rF(v+x)  -rl. 

A  ■  kx*  -  7  x(e  '«  t  ~  P 


F(v+x) 

|  (D(y) 


v)e’rydy 


The  integrand  is  bounded  as  follows: 


S3 


0  «  y  «  V(v  ♦  x)  , 
0(0)  <  D(y)  *  v  ♦  x  , 
0(0)  -  v  «  !>(y)  -  v  <  x  . 


Substituting  the*  lower  bound,  gives: 


«  ki*  -  e  -  .-rI)  *  *  (m(o)  -  v)(.'r''(v**,  -  I)  . 


As  in  Case  111,  the*  right  side  of  this  inequality  Is  negative  whenever 


(1  -  a)kx-  <  ?  (0(0)  -  v)(l  -  e"rP(v4x))  , 


or 


(10a) 


u  -  •>(*)'  <  a  -  . 


Hence,  as  in  Case  III: 


(1  -  a)w1_axa  <  (D (0)  -  v)(l  -  e"rF(v+x))  ■>  T(x)  <  T(0) 


Substituting  x  ,  the  upper  bound,  wc  have: 


4  >  k;*  -  h  ;<e-rr(vt,!)  -  *-ri)  -  £  ;<i  -  «-rF<vt,)) , 


» »;*[i  -  (if- « - .-*>] . 


Again  as  in  Case  III,  the  right  side  of  this  inequality  is  positive  whenever: 


erI  < 


-  r 1 


or 
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(10b) 


hence  s 


'I  <  i  tn 

■  r 


1 


-  r 1 


->  T(x)  >  T(0)  , 


the  earn*  at  we  found  for  Cote  III.  If  neither  (10a)  nor  (10b)  holds,  a 

* 

direct  comparison  of  T(x)  with  T(0)  mutt  be  made. 

Cane  IV:  Summary 

(a)  x  c  (O.w)  . 

t*(x)  -  I  and  T(x)  it  concave  increasing  for  all 
x  c  (0,w]  . 

(b)  x  c  (w.D(I)  -  v]  . 


itHrr^ 

V  ■  a\D(I)  - 


— 1  -O  T’(x)  >  0  V  x  e  (w,D(I)  -  v] 
,l-a|  * 


and,  hence,  x  ■  0  . 

But  when  this  inequality  does  not  hold,  there  may  exist 

a 

relative  minima  for  x  e  (D(0)  -  v,D(I)  -  v)  ,  and  x  yields  the 
smallest  value  of  T(x)  in  this  interval.  But  x  may  not  be 
the  optimal  value,  unless  T(x)  <  T(0)  .  And: 


x  *  x  if  (1 


x  -  0  if  I  <  r 


<  (1  -  e 


-  ■>© 

•l‘M 


-rF(vfx) 


) 


and  if  neither  of  these  inequalities  hold,  then  T(x)  must  be 
compared  directly  with  T(0)  . 


MfrtHAKV  OF  At.l  CASKS  M>K  omHAI.  k 
li  w  >  {0(1}  -  v)  . 


*  „  _  * 

t  (x)  »  1  ,  x  •  0  . 


Casc_  IJj _ |i(0)  -  v  <  w  <  P(jl)  -  v 


/  *  „ 
x  ■  0 1 


t  (x)  -  I 


Let  T(x  )  be  the  binal'est  value  of  T(x)  where  x  satisfies 
o  o 


1-a  a 


akx 


a-1  -rF(v+w"  x  )  p  .  -rF(v+x)  -rl. 
e  ■  *•  (e  -  e  )  . 


Let : 


U(xo,v,r,w,a) 


L(xo,v,r,w,a) 


Then: 


*  \ 

x  -  x  I 

°  ( 

>  if  I  >  U(x  ,v,r,w,a) 

t*(xo)  -  f(v  +  w1_V)j 


and  if  L(xQ,v,r ,w,a)  <  I  <  U(yQ  ,v,r ,w ,a) 


■> 


a  direct  comparison  is  required. 


Bacall: 


(8) 

(8a) 


(9) 


(9a) 


(U) 


If  (11)  and  (8a)  hold,  then  x*  -  0  because  T(x)  is  nondecreasing 
throughout  the  entire  interval. 

If  (11)  and  (9)  hold  there  are  no  relative  minima  in  (z,D(I)  -  v)  but 
there  may  be  relative  minima  in  (D(0)  -  v,z]  ;  and,  if  there  are: 


and  T(x)  must  be  compared  with  T(0)  if  neither  of  these  tests  apply, 


in  which  case, 
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* 

x  •  x 


t  (x)  •  0 


if  T(x)  <  T(0)  , 


(x  -  0}  if  T(x)  >  T(0)  . 

If  (11)  does  not  hold  but  (8)  or  (8s)  holds,  there  ere  no  reletlve 
minims  in  (v,z]  but  there  mey  be  reletive  minima  in  (t,D(l)-v);  end  if 
there  ere: 


* 

x  ■  x 


k<*.>  -  * ♦ 


if  1  >  U(xo,v,r,w,s)  , 


(x  -  0}  if  I  <  L(xo,v,r,v,a)  , 


but  if  L(xQ,v,r ,w,s)  <  I  <  U(xo,v,r,w,a)  ,  a  direct  comparison  must  be 


made  and 


* 

x  ■  x 


[t*(xo)  -  f(v  +  v1_V)j 


If  T(xq)  <  T(0)  , 


{x  -  0}  if  T(x  )  >  T(0)  . 

o  ■ 


If  (9)  holds  but  not  (11)  there  may  be  relative  minima  in  both  Intervals, 
(D(0)  -  v,z]  and  (z,D(I)  -  v]  and  these  must  be  compared  vith  T(0)  .  And, 


(  x*  .  X  1 

I 

-  if  ' 

^T(x0)  <  T(0)  and| 

(t*(xo)  -  f(v  +  v1’8**) 

J(xo)  <  T(x)  j 

V 
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X  ■  X 

t"(x)  ■  o( 


if 


and 


(x  ■  0)  if 


T{x)  <  T(0)  and 
T(x)  <  T(x) 


T(0)  <,  T(x  )  and 
o 

A 

T(0)  <  T(x) 


Case  IV:  W  <  D(0)-v  and  z  £  D(l)-v. 


Recall: 


(12) 


1  .  (  “ 
<  r  ’  /n  I  1  - 


trcITT^) 


If  (12)  holds,  there  are  no  relative  minima  in  the  interval 
(w,  D(l)-v]  and  x  «0.  Otherwise,  there  may  be  relative 
minima  in  the  open  interval,  (v,D(l)-v).  But,  x,  the 
value  of  x,  for  which  T(x)  is  a  minimum  in  this  interval, 
i3  not  the  optimal  value  unless  T(x)  <  T(0). 

Then: 

x*  =  x  if  (l-a)(^)  <  (l  _e-rF(v  +  x)) 
x*  =  0  if  I  <  I  In 

a  r 


and  if  neither  of  these  tests  are  met,  a  direct  comparison 
of  T(x)  with  T( 0)  is  required. 
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Remarks  on  uniqueness  of  relative  minima  x  and  x 


2 


Zn  our  analysis  of  Cases  II ,  III  and  IV  ve  found  that  whenever 
relative  minima  exist  on  the  interior  of  x  c  (w,D(l)  -  v),  they  are 
for  values  of  x  which  satisfy  T'  (x)  ■  0,  or: 


T'(x)  -  xkx*-Vrt‘(x)  -  .  .-■’I)  .  o 

r 

Accordingly,  denote  S(x)  for  x  c  (v,D(l)  -  v)  as  follows: 

S(x)  -  Hx*-1  ,-rt’(x)  -  I  (e-rP(v+x)  -  e-rI), 

r 

where:  . 

IF(v  ♦  v1_axa)  for  x  e  (w,D(l)  -  v]  for  Case  II, 

F(v  ♦  w1_axa)  for  x  c  (z,D(l)  -  v]  for  Case  III, 

0  for  x  c  (D(0)  -  v,z]  for  Case  III, 

0  for  x  e  (D(0)  -  v,l]  for  Case  IV 


Denote : 


H(x)  -  akxa-1e~rt"(x), 

G(x)  •  »  (e-rF(v+x)  .  e-rI)  . 
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Them 

S(x)  -  H(x)  -  G(x)  , 

S '  (x)  -  H'<x)  -  C* (x)  , 
and 

S"(x)  -  H"(x)  -  G"(x)  . 

Then: 

H' (x)  ■  .(.  -  l)kx*'Vrt*(x)  -  r.kx-le-rt*(x)(^  t*(x))  . 

And 

HM(x)  -  a(a  -  1)  (a  -  2)kxa“Vrt  (x)  +  (1  -  a)rakxa"2e"rt  (x)[^  t*(x)] 

+  rW-V1'*00^  t*(x)]2  -  r.kx-1«-rt‘<x>^  .*«]  , 

2 

q  *  d  * 

which  is  positive  whenever  — j  c  (x)  is  negative  or  zero  because  ^  t  (x) 

dx* 

it 

is  nonnegative.  (H'’(x)  >  0  for  t  (x)  ■  0)  .  Now: 

G'M  -  pe-rF(v+x)  ±  F(v  +  *)  »  0  . 

And: 

G-W  -  ,.-*<«*>/-  F(y  +  x)]2  +  JL  F(v  +  x)j  , 

d2 

which  is  nonpositive  whenever  — r  F(v  +  x)  is  nonpositive,  which  in  turn 

dx^ 


implies  -G'^x)  is  positive. 
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Hence  S"(x)  la  poaltivo  and  S(x)  la  atrlctly  convex  whenever  both 
d2  d2  * 

— x  F(v  +  x)  and  — x  t  (x)  are  nonpoaltlve.  One  clasa  of  demand  functions 
dxz  dx 

for  which  thaae  condition*  hold  la  D(y)  convex  becauae  F(x)  la  concave. 
For  auch  a  claaa,  thare  are  at  moot  two  aeroa  of  S(x)  . 

For  example,  auppoae  v  -  0  ,  D(y)  *  ey  and  w  ■  1  .  Then  Caae  II 
applies  and  t*(x)  -  F(v  +  w2_axa)  .  Restrict  x  c  (w,D(I)]  •*>  1  <  x  <  e*  . 
Then 

in  D(y)  -  y 

F(D(y))  -  F(ey)  -  in  D(y)  -  y  , 

F(v  +  x)  ■  in  (v  +  x)  *  in  x  , 

£  Flv  +  x>  •  I  • 

t*(x)  -  F(xa)  »  a  in  x  , 

d  *.  .  a  d2  *,  .  -2 

to  c  (x)  ‘  x  ■  ~2  1  (x)  -  -*x 
dx 


— x  F(v  +  x)  -  -x  <  0  . 
dx 


Further: 


H''(x)  ■  ka(a  -  1) (a  -  2)xa  ^e  rt  ^  +  (1  -  a)rakxa  2e  xt 


2  a-1  -rt  (x)/a\  ,  a-1 

+  r  ax  e  I  — 1  +  rax  e 


-l--rt  <XW2>  . 


Hence  H'  (x)  >0  .  Now: 


G' '  (x)  -  pe 


-rinx 


(-&  ■  ■") 


<  0  . 
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» 


Hence  -C-"(x)  >  0  f  S"(x)  -  H"<x)  -  <•"<»  '  U  and  S(x)  is  convex 
with  at  most  two  zeros. 

For  Case  II  and  IV  when  t*(x)  -  0  for  all  x  t  (w,z)  : 

H(x)  |  c  akxa  1  >  0  , 

|xe(w,z) 

and  H(x)  is  convex,  decreasing.  Hence  these  remarks  apply  to  all  cases 
except  Case  I,  for  x  e  (w,l)(l)  -  vj  . 


An  Algorithm  for  Finding  an  Optimal  Expansion  Policy 


1.  Define  whether  Case  I,  II,  III,  or  IV  applies  as  follows: 
Compute  w  -(D(l)-v).  If  non-negative,  Case  I  applies. 
Otherwise,  compute  w  -(D(O)-v).  If  non-negative.  Case  II 
applies.  Otherwise,  compute  D(l)-v-z.  If  non-negative, 

Case  III  applies;  otherwise  Case  IV  applies.  If  Case  I 

#  . 

applies,  x  =  0.  For  Cases  II,  III,  IV,  go  to  stops  2,3,4 

respectively. 

2.  (Case  II).  Compute: 

a/  w 
\D(l)-v 

If  non-negative  x  =  0  because  T(x)  has  no  relative  minima 
exist  in  the  interior  of  (0,D(l)-v). 

Otherwise  go  to  step  5- 


1-a 


-rl  -rF(v+w) 
+  e  -e 


3*  Case  III.  Compute  A  and  A  where: 

1  2 


A  =  a 
1 


/  .  ya 

(d(i)-v  j 


-rF(v+z)  -rl 
e  -  e _ 

-rF(v+w1_a(D(l)-v)a) 


and 


A  =  a  -  /z{"a  (l-e"rI) 
2 


(i) 
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* 


Remark : 

A^  >  0  ■>  inequality  ( 11 )  holds,  and 

A  >  0  =>  inequality  (9)  holds  =>  (10)  does  not  hold. 

2  “ 

*  ,  . 

If  A  >  0  and  A  >  0,  x  B  0  because  T(x)  has  no  relative 

1=2  = 

minima  in  (0,  D(l)-v]. 

If  A^  <  0  and  A^  >  0,  there  may  be  relative  minime ,  but  only  in 

(z,  D(l)-v);  hence  go  to  step  5. 

If  A  >•  0  and  A  <0,  there  may  be  relative  minima,  but  only  in 

1  “  2 

(D(O)-v,  z];  hence,  go  to  step  6. 

If  A  <  0  and  A  <0,  there  may  be  relative  minima  in  both  intervals 

1  2 

(D(0)-v,z]  and  (z,  D(lO-v);  hence,  go  to  step  7. 


1*.  Case  IV.  Compute 


If-  non-negative,  x  =  0  because  T(x)  has  no  local  minima  in  the 
interior  of  (0,D(l)-v).  Otherwise  go  to  step  6. 

5.  (Cases  II  and  III).  Compute  all  solutions  to: 
a(*)a  l(  e-rFU^-Vlj*  (e-rIr(T+x)-  <frI) 

where:  x  e(v,D(l)  -  v)  for  Case  II  and  x  c  (z,  D(l)-v)  for  Case  III. 
Newton's  Method  or  some  other  appropriate  method  of  numerical  analysis 
may  be  employed.  If  D(y)  is  convex,  there  are  at  most  two  solutions. 
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Denoting  xq  as  the  "best"  solution,  compute 
I  -  L(xo»  v,  r,  v,  a)  where 


« 

If  the  result  is  positive,  x  =  xq  .  Otherwise,  compute 
I  -  U(xo,v,r,v,a) 

where : 

U(x  ,v,r,v,a)  =  F(v  +  w1-axa)  +  —  In 
o  or 

If  this  result  is  negative  or  zero,  x  =0.  Otherwise, 

£ 

compute  T(xq)  -  T(0).  If  this  result  is  negative,  x  =  0; 

* 

otherwise  x  =0. 

6.  (Cases  III  and  IV).  Compute  all  solutions  to 


(x\  1-a  -rF(  v+x)  -rl 

%)  e  "  e 

where  x  e(D(0)-v,  z]  for  Case  III  and  x  e  (D(0)-v,  D(l)-v) 
for  Case  IV. 


If  D(y)  is  convex,  there  are  at  most  two  solutions.  Denoting 

A 

x  as  the  best  solution,  compute 

(l-a)  (f)a  -  (l-e-rr(v+x)). 

If  the  result  is  negative,  x  =  x.  Otherwise,  compute 


v» 
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» 


If  the  result  is  non-negative,  x  =0.  Otherwise,  compute 

* 

T(x)  -  T(0).  If  the  result  is  negative,  x  =  x,  otherwise 

. « 

x  =  0. 

7.  (Case  III:  both  x  and  xq  may  exist.)  Compute  all  solutions 
to 

e-rF(v+w1_axa)  =  (e-rF(v+x)  _  e-rl)# 

where  x  e  (z,  D(l)-v),  and  xq  is  the  best  solution.  Compute 
all  solutions  to 

a/^a_1  .  (e-rF(v«)  .  ,-rIj  _ 

where  x  e  (D(0)-v,  z]  ,  and  x  as  the  best  solution.  Compute 

A 

T(x  )  -  T(0 ) ,  and  T(x)  -  T(0). 
o 

T(0)  and 

A 

T(x) 

T(0)  and 
T(x0) 


Then: 


* 

x  =  x. 


* 

x'  -  X 


if 


if 


T(x  )  < 
o 

T(x  )  < 
0 

T(x)  < 
T(x)  < 


x  =  0  otherwise, 
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CHAPTER  IV 

AN  ALGORITHM  FOR  SEQUENTIAL  CAPACITY  EXPANSIONS 

Introduction 

In  this  chapter  we  consider  a  sequence  of  capacity  expansions  over  time. 
The  algorithm  of  Chapter  III  for  a  single  expansion  is  too  cumbersome  to  use 
for  sequential  expansions.  Accordingly,  we  present  an  algorithm  for  digital 
computers  in  which  demand  and  capacity  take  on  discrete  values. 

Our  approach  is  to  treat  an  n-expansion  problem  of  2n  decision 
variables  as  an  n-stage  decision  process.  Two  types  of  problems  are 
considered.  The  first,  Type  I,  does  not  permit  a  shortage  of  capacity  after 
the  last  expansion  while  a  Type  II  problem  does.  We  show  that  solving  a 
Type  II  problem  requires  that  a  Type  I  problem  also  be  solved. 

Let  Xj  denote  the  size  of  the  Jth  capacity  expansion,  t^  the  time 
at  which  it  occurs  and  v^  the  capacity  before  it  occurs.  Restrict  the 
times  of  capacity  expansions  as  follows: 

0  <  t,  <  t_  <  ...  <  t  <1. 

■  1  "  2  -  ■  n  « 


where  I  is  the  length  of  the  planning  interval  and  t  e  [0,1]  . 

Suppose  we  restrict  x  such  that  x  ■  D(I)  -  v  and  allow  (D(I)  -  v) 
to  vary.  This  forces  the  last  term  of  T(x)  to  vanish,  and 
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I 

(13a)  p  J  max  (D(y)  -  v  -  x,0)e  rydy  =  0  . 

t 

Let  C(«)  denote  the  single  expansion  cost  function  under  this  restriction. 
To  avoid  ambiguity,  let  the  last  expansion  for  a  j-stage  problem  satisfy: 

-  *><*>  -  *}  ■ 


where 


and 


D(i)  e  [D(0),D(I>] 
Vj  e  [ vQ  ,D(I)  ]  . 


Note t we  are  also  defining 
capacity  at  time  zero. 


to  be  a  parameter,  and  vq  to  be  the  initial 


Remark: 

A 

Note  that  by  our  restriction  on  x  ,  t  (x)  <  F(v  +  x)  .  Hence,  the 
configuration  implied  by  our  definition  of  Type  I  problems  is  shown  below 
(for  <  D(0))  . 


FIGURE  7 


69 


In  particular,  we  can  write  without  loss  of  generality  that 

&  A 

0  <  t^x^)  <  F(v^  +  x^)  <  t2(x2)  <  F(v2  +  x2)  <  ...  <  I  . 

This  will  permit  us  to  simplify  the  cost  functions.  Note  also  that  we  can 
define  the  jth  stage  unambiguously  as  beginning  at 

Vj  *  F^V1  +  xi  +  •**  +  xj-l^  an<*  encIin8  Vj+1  *  FCv^  +  x^  +  ...  +  Xj)  . 

Restricting  x  ,  so  that  (13a)  holds,  the  cost  function  for  the  Type  I 
single  stage  problem  is 

*  *1^*1^ 

-rt.(x  )  /• 

C1^X1,VP  "  kxle  +  P  I  max  (D(y)  -  vlt0)e“rydy  , 

0 

where : 

-  (D(i)  -  Vj)  , 
and 


t*(xi) 


2  2 


if  x^  E 


'  A 

F(vl  +  3xl)  otherwise 


(w,z)  and  w  +  <  D(0)\  . 


Note  that  t2(x2)  <  I  with  equality  always  holding  for  a  one-stage  problem. 

The  cost  function  for  the  first  two  stages  of  a  Type  1  problem  is 
given  by  the  sum: 


* 
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Cl(xl’vl)  +  C2(X2,V2)  "  ^ 


,  j 

o 


1 


max  (D(y)  -  v1,0)e'ry<ly 


-rt0(x,) 

+  kx^e  +  P 

?«*!> 


max  (D(y)  - 


v2,0)e-rydy 


where 


:1  +  X2  *  “  V1  ’ 


t2(x2)  = 


v2 

-  vl  + 

ft*3(x3) 

if  x2 

{*>!> 

/  / 

if  x2 
1-a  a\ 

Y  \v2  + 

w  x2) 

and  tt(x^)  -  I  if  this  be  only  a  two-stage  problem.  This  can  be  simplified 

of  definition.  Write  C^.v^  in  an  equivalent  form: 


'3'  3 

to  obtain  symmetry 


:i(x1.,1)  -  kx^e"rt‘<'tl)  +  P  JVl>  -  CDCy)  -Vj,0)e_tydy 


+  p 


j  max  (D(y)  -v^,0)e  rydy 


F(vx) 


Examine  the  three  cases  for  t^(x^)  • 


Case  1: 


t*(xx)  -  0  «>  F(vx)  -  0  . 


Then  C^(x^,v^)  may  be  written  in  equivalent  form  as 

*  * 

-rt.(x.)  r 

C1^X1,V1^  *  kxle  +  p  I  0>(y)  “vl»°^e~r 

F(vx) 

Case  2 : 

t^x^)  "  F(vl  +  wl  Sxl)  >  F^vi^  * 

Now,  for  the  integrand  of  the  first  integral: 

0  <  y  <  , 

D(0)  <  D(y)  <  vx, 
or  D(y)  -  vx  <  0  . 

Hence  max  (D(y)  -  v^.O)  *  0  . 

And,  for  the  second  integrand 

F(vj)  <  y  <  f(vx  +  w1_ax®)  > 

vi  m  °(y)  > 
or 

0  <  D(y)  -  v1  . 

Hence  max  (D(y)  -  v^,0)  *  D(y)  -  . 

Again  C^Cxj^jV^)  may  be  written  as  above  in  Case  1. 


Case  3: 


Again  as  before,  for  y  e  [0,F(v^)j  ,  max  (D(y)  -  v^,0)  ■  0  ,  and  for 
y  t  Jf^)  ,t*(x2)J  ,  max  (D(y)  -  v^O)  -  D(y)  -  . 

Hence  without  loss  of  generality  we  can  write 

*  **<*!> 

-rt.^)  r 

(14)  »  kxae  +  p  I  (D(y)  -  v^e  ydy 

F(v1> 

where  F(v,)  *  0  whenever  v.  <  D(0)  . 

1  Is* 

Remark : 

A  A 

Case  3,  t^(x^)  *  t2(x2)  is  a  degenerate  case  but  presents  no  real 

difficulty.  Whenever  two-stages  "collapse"  into  one  stage,  the  cost  function 
for  a  single  stage  is  less  than  for  two  stages,  for  expansions  at  the  same 
point  in  time.  As  an  example,  for  v2  “  vi  +  xi  and  ^  5 

CCx^v^  +  C(x2,v1  +  x2)  >  C(x1  +  x2,v1) 

because 

(kx*  +  kx2j  >  k(x.^  +  x2)a 

for  0  <  a  <  1  . 

( 

Now,  consider  the  second  stage  cost  function,  and,  as  before  examine 
A  A 

the  last  term.  If  t2(x2)  *  ti^xl^  IflSt  terra  Vilnis^es>  If 

A  A 

t2(x2)  *  t3(x3)  the  case  is  degenerate  and  three  stages  have  collapsed  into 

*  /  1-a  a\ 

two  stages.  However,  if  t2(x2)  *  F(v2  +  w  x2  )  then  the  arguements  similar 
to  those  leading  to  (14)  can  be  made.  Write  the  integral  as  follows: 
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t* (x2) 


FfvJ 

4. 


J  max  (D(y)  -  Vj.OJa-'J'dy  •  p  J  max  (D(y)  -  Vj .OJ.'^dy 


ti(x1) 


tl<xl) 


t2^x2> 


+  P  J  max  (D(y)  -  v?,0)e"rydy  , 


F(v2) 


where  v2  -  ^  +  xx  and  P(v2)  -  F(Vl  +  x^  .  Then,  the  first  integral 
always  vanishes  because 

tl<xi)  <  y  <  F(vx  +  xA) 


and 


D(y)  <  v1  +  Xx  , 


or 


D(y)  -  vx  -  xx  <  0  , 


Hence 


max  (D(y)  -  v2>0)  -  0  . 


Also,  for  the  second  integral. 


F(v2)  -  F(Vl  +  xx)  <  y  <  t*(x2)  , 


or 


V1  +  *i  <  D(y)  . 


♦ 


Hence : 


max  (D(y)  -  v2,0)  ■  D(y)  -  v2 


Therefore,  the  cost  function  for  the  second  stage  can  be  written 


*t  V  t2(x2) 

-rt  (xj  r 

c2Cx2*v2^  "  kx2e  +  P  J  (D(y)  -  v2)e  ^dy  , 

F(v2) 


9* 

where  t2(x2)  is  defined  as  before.  In  general,  then,  we  can  write  for  the 
Jth  stage. 


(15)  C 


*  *!<*4> 

a  •rti<xj)  /  J 

"  kXje  +P  J 

F(v.) 


(D(y)  -  Vj)e-rydy  , 


where:  Vj  -  Vj  +  xx  +  . . .  +  . 


t*(x  )  »  0  i  and 
0  o 


(16)  tj(Xj)  “ 


if  xj  c  1° I 

t^_1(xj_1)  if  Xj  c  (w,z)  and  w  +  Vj  <  D(0)5  • 


F^Vj  +  w1  otherwise 


Computational  Procedure,  Type  I  Problems 


1.  Divide  the  capacity  and  demand  interval  [0,D(I)}  into  M  subintervals. 
For  illustrative  purposes  only,  assume  these  subintervals  are  of  uniform 
size,  say  A  >  0  .  Consider  v^  ,  the  initial  capacity  to  be  a  parameter 
instead  of  a  fixed  value,  with  v^  e  ( vq , D ( X ) ]  where  vq  and  D(I) 
are  fixed  values.  Restrict  positive  values  of  ,Xj  ,  D(*)  ,  Vj  to 
multiples  of  A  .  Hence  ,  for  a  J-stage  problem, 
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v  +  x.  B  D(i)  e  {0,A,2A,  ...,MA) 
J  J 


v  c  {0,A,2A,  .  ..,  MA} 
J 


(Note:  except  for  linear  demand,  F( ‘ )  is  not  restricted  to  multiples 
of  some  value.) 

2.  Choose  a  value  for  D(i)  and  a  value  for  v^. 

3.  Compute  t*(x^)  from  (13). 

I*.  Compute  C1(x1>v1)  from  (ll*). 

5.  Repeat  Steps  2, 3, and  It  for  all  D(i)  c  [v0,D(l)]  and  all  v^c[votD( I ) ] . 

6.  Store  the  results  in  Table  I.  The  subscript  J  is  used  because  Table  I 
is  used  for  any  succeeding  single  stage  cost  function. 


Cj(  D(i ) ,Vj  )  J*l*2 . n 


D(i )  v0  v0  +  A  vc+2A  .  D(I)-2A  D(l)-A  D(l) 


v0+2A 


.1)  -  2A 


,D  -  A 


TABLE  I 

(Note:  diagonal  entries  are  zero  and  all  entries  below  the  diagonal 
are  also  zero). 
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7.  Now,  consider  a  two  stage  problem.  The  first  stage  end  conditions 
are  the  second  stage  initial  conditions.  Therefore,  demand  equals 
capacity  at  the  beginning  of  the  second  stage,  except  when  x^  =  0  or 
x^  =  D(l)-v0,  in  which  case  there  is  no  second  stage.  Hence  both 
first  and  second  stage  cost  functions  are  found  in  Table  I.  For 
example,  suppose  x^=  2A>w  and  =  v0  for  the  first  stage,  and. 

x0  =  D(I)  -  2A  -  vQ  v  =  v0  +  2A  for  the  second  stage.  Then,  the 
entry  in  Table  I  for  the  second  stage  cost  function,  corresponds  to 
column  D(l)  and  row  (v0+  2A).  For  this  example,  the  asterisk  in  Table  I 
is  for  the  first  stage  and  the  double  asterisk  for  the  second  stage. 

8.  Restrict  the  sum  of  x^  and  x^  such  that 

X1  +  x2  =  [  D(!)  -  v°  ]. 

Let  =  v0  .  Then,  since 

v?  =  vQ  + 

Table  I  can  be  used  to  compute  by  enumeration: 


C*(D( I ) )  =  minimum  {  C0(D(l),v0)  +  Ctf (d( I )  -  v_)  }, 
2  v2c[v0,D(I)]  221  2 


where  C*( D( I ) )  is  the  optimal  solution  for  a  one  stage  problem. 

This  yields  the  optimal  policy  for  a  capacity  expansion  of  size 

D(l)-v0  done  in  n=2  stages.  But  if  expansion  in  two  stages  is  not 

« 

as  good  as  expansion  in  one  stage,  x^  =  0.  Hence,  this  also  yields 

* 

the  best  value  for  n  ,  the  optimal  number  of  stages. 

9.  Now,  consider  a  three  stage  problem,  with  a  third  stage  initial 
capacity  of  v^,  where: 


V3 


ro  +  xx  +  x2 
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For  a  given  value  of  v^  «  the  optimal  total  cost  function  is  given 
by  the  sum  of  the  third  stage  cost  function  plus  an  optimal  two 
stage  cost  function,  i.e., 

C3(D(I),v3)  +  C*(D(I)  -v3). 

We  wish  to  minimize  this  sum  by  choosing  the  best  value  for  v^. 

# 

Hence,  compute  by  enumeration,  Cg(D(l)  -  v  )  for  all  v^  e  [v0,  D( I ) ] 
and  store  in  Table  II. 


V3 

C*(D(I)-v3) 

0 

A 

• 

• 

• 

D(I) 

TABLE  II 

10.  Now,  use  Table  I  and  II  to  compute  by  enumeration: 

C,(D(l))  =  minimum  {  C_(D(l),v  )  +  C  (D(l)  -v  ) }  . 

3  v3e[v0iD(l>]  332  3 

This  yields  the  optimal  policy,  including  n  ,  for  an  n=3  stage 

problem  and  if  expansion  in  three  stages  is  not  as'  good  as  two  stages , 

*  . 

=  0.  Uf  expansion  in  three  stages  is  not  as  good  as  one  stage, 

«  « 

x3  =  x2  =  0. } 
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11.  Repeat  step  9  and  10,  adding  a  fourth  atage  and  construct  a  Table  III 
of  values  for  C*(D(I>-v,,),  using  Tables  I  and  II.  Discard  Table  II 

if  desired. 

when  the  previous  table  is  repeated.  If  this 


12.  Terminate  the  process 


happens 


for  n=k+l ,  then  n  =  k. 


Remark : 


We  can  now 


write  the  general  recursion  relation. 


(20)  C*(D( I ) )  =  minimum  j  ~Vn  +  (^(Dd)  .VR)  J 

n  v  e[v0,D(l)J  ( 

nc 


Type  II  Problems 

Consider  now  a  type  II  problem  in  which  the  shortage  after  the 
last  expansion  is  permitted.  Denote  the  shortage  after  the  nth 
capacity  expansion  as: 

I 


s(vn)  =  p  J  -  \  )e“,dy- 


F(vn) 


This  function  is  the  no  expansion  case  whenever  vn  =  v„  . 

Denote  T2(  • )  as  a  two  stage  Type  II  cost  function.  For  a  given 
value  of  v2,  the  two  stage  cost  function  msy  be  written  as  the  sum  of 
the  shortage  cost  plus  an  optimal  n  "2  stage  Type  I  cost  function: 


T2(v2)  =  C*(  v2)  +  S(v2) 
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We  wish  to  minimize  this  by  choosing  the  best  value  for  Vg.  Hence,  we 
define 


*  (  • 

T0(D(l) )  =  minimum  VC  (v_)  +  S(v  ) 

2  v2  e  [v0,D(l)] I  2  2  2 

as  the  optimal  cost  function  for  a  two  stage,  type  II  capacity  expansion 
problem.  Similarly,  for  an  n-stage  problem  we  have: 

T  (d(I))  =  minimum  I  C  (v  )  +  S(v  ) 

"  vn  E  [v0,E(l))|  n  n 

The  Type  II  computational  procedure  parallels  the  type  I  procedure 
after  a  table  of  values  for  S^)  has  been  computed  for  all  v^  e  [vOJD(I)]. 

(Note  that  the  no  expansion  policy  is  not  found  by  type  I  computations 
but  it  is  found  by  type  II  computations.) 
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CHAPTER  V 


A  STATIONARY  POLICY  FOR  A  LINEAR  DEMAND  FUNCTION 


Introduction 


In  this  chapter  we  use  our  previous  results  for  a  finite  sequence  of 
capacity  expansions  and  argue  an  asymptotic  stationary  model  as  a  limiting 
case  for  a  linear  demand  function,  with  a  planning  interval  of  infinite 
length  and  an  infinite  number  of  stages.  Manne  [9]  who  first  investigated 
this  problem  uses  the  notion  of  "regeneration  points"  to  derive  a  cost  function. 

Consider  a  sequential  capacity  expansion  problem  in  which  the  demand 
function  is  linear  and  the  initial  demand  as  well  as  initial  capacity  are 
zero.  Let  the  rate  of  growth  of  demand  be  "b"  .  Then 


D(y)  *  by  and  D(0)  -  v  =  0 


F(x)  -  ^  x  . 


From  Equation  (16),  Chapter  IV  we  know  that  for  vq  *  D(0)  ■*  0  ,  the  optimal 
time  of  expansion  for  the  nth  stage,  n  «=  1,2,  ...,  N  ,  reduces  to  two  cases: 


t  <X>  E 

n  n 


t  (x  , if  0  <  x  <  w 
n+1  n+1  -  n  = 


f(v  +  w*  V  otherwise 
.  V  n  n/ 


where 


n-1 

v  ■  v,  +  y  x. 
n  1  i 
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Also  from  Chapter  IV  we  know  the  configuration  of  the  optimal  expansion 

policy  is  F(x,  +x„+...+x  , )  <  t  (x  )  <  F(x.  +  ...  +  x  )  because 

1  l  n-1  ■  n  n  ■  1  n 

■  0  .  Hence  we  may  write  the  cost  function  for  an  n-stage  Type  I  problem 

as  follows.  Let:  x  ■  t  B  0  and  t  ■  I  .  For  brevity  denote: 

o  o  n+1 


C 

n 


cl(Xi,Vi)  +  C2(x2,v2)  + 


+  C(x,v) 
n  n  n 


Then: 


C 

n 


(by  -  v^e  rydy  . 


Define  a  stationary  policy  as  one  in  which  all  capacity  expansions  are  of 
equal  magnitude.  This  is  equivalent  to  specifying  that  the  time  intervals 
between  expansions  are  equal  for  an  optimal  policy.  To  see  this  consider 
that  the  time  interval  between  the  nth  and  n  +  1th  expansions,  is: 


And: 


x 


1 


Hence 
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*  *  1-a  a.  ,v  .  1-a  a> 

t  ,,  -  t  B  F(nx  +  w  x  )  -  F((n  -  l)x  +  w  x  ) 
n+1  n 

1  /  ,  1“3  a.  1  §  j  « v  ,  1“ a  a. 

•»  r-(nx  +  w  x  )  -  --( (n  -  l)x  +  w  x  ) 
b  b 

1,  ,  1-a  a. 

b 

and  the  time  intervals  between  successive  expansions  is  constant.  In  particular 

*,  s  1  1-a  a 

cl(x)  *  £  W  X 

t*(x)  »  F (x)  +  t*(x) 

and 

t*(x)  =  (i  -  l)F(x)  +  t*(x)  ;  i  =  1,2,  ....  n  . 

The  asymptotic  cost  function  as  n  -*•  °°  and  I  “  is  denoted 
and  using  the  preceding  we  may  write  it  as  follows. 


00  _rt  00  (  iF(x)+t1  | 

C  «  l  kxae  1  +  £  <p  f  (by  -  ix)e  rydy >  . 

1-1  1=0  (  iF(x)  j 


Note : 


~rti+l  I  -r(iF(x^+ti 

e  ■  le 


’]  -  [e’^W]1. 


-rt. 


hence 


£  kx“e  *  =  \kx“e 
i=l 


8  rtl  -  (kx8e'rtl){l  +  e-rF(x)  +  [e-rF(x)]2  +  ...  } 


-  U^)  i  [.-"W]1 

i-0 


-rt. 
kxae  1 


1  -  e 


-rF(x) 
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-rF(x) 

became  |e  <1  for  x  >  0.  Next  examine  the  second  term,  the 


integral,  in  . 

iP(x)+t. 
P 


iP(x)+t1 

f 

iF(x) 


bye-rydy  -  pix 


iF(x)+t. 

j  «_rydy 

iF(x) 


f 

=  e  riF^x^  | [b  +  briF(x)  -rix]  +  e  (rix  -  b  -  briF(x)  -  brt^) j 


To  compare  this  result  directly  with  the  asymptotic  stationary  model  of 
Manne  [9l>  we  normalize  the  time  scale,  letting  b  =  1  and,  therefore, 
F(x)  =  x.  Then  the  above  expression  becomes 


h 


P _  -rix  I  ,  -rt-i  /, 

— 2“  e  |  1-e  1  (1  +  rt 

r 


.>]• 


Then 


C.  *  ’f  '-rx1 '  *  4  (  ^  rtl  (1  *  rtl)£fe-rX]  1 

1-e  r  \  'i=o 


-rx 


Again,  noting  that  |  e  |  <  1  for  x  >  0  ,  we  can  write 

iA-Y-1 


i=0 


1-e 


-rx 


Hence 


(^]["0e'rtl  +  ?(u-e’rtl(1  +  r^)]  ' 


This  result  agrees  with  Manne  [9].  It  can  be  simplified  by  noting  that 


* 
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t*(x) 


rk  a 
—  x 
P 


or 


£ 

r 


tjCx) 


Substituting  this  into  the  above  yields  a  cost  function  in  one  variable,  x  : 


C  (x) 
00 


This  is  the  same  result  as  that  found  by  Erlenkotter  in  Chapter  11  of 
Manne  [lO].  There  the  derivation  was  based  on  the  notion  of  "regeneration 
points"  rather  than  directly  arguing  the  case  as  an  asymptotic  result  of  a 
sequential  capacity  expansion  model. 

In  the  chapter  which  follows,  we  will  use  the  algorithm  of  Chapter  IV 
to  compare  a  finite  planning  interval  policy  with  a  stationary  policy. 


Remarks ; 

There  are  distinations  between  assuming  an  infinite  planning  horizon  and 
a  finite  planning  horizon.  One  reason  Manne  suggests  an  infinite  horizon  is 
that  he  doubts  that  geometric  growth  rates  which  characterize  some  industries 
today  "can  be  extrapolated  to  the  indefinite  future."  Another  reason  he 
gives  is  that  his  cost  function  is  mathematically  tractable.  Finally,  his 
model  sheds  light  on  the  question  of  optimal  policies  for  standard  sized 
plants.  However,  the  chances  of  correctly  forecasting  a  demand  function 
far  into  the  future  are  very  small  indeed,  hence  a  policy  based  on  such 
deterministic  forecasts  is  not  truly  optimal.  On  the  other  hand,  a  policy 
based  solely  on  short  term  forecasts  may  not  be  truly  optimal  either. 

Although  it  may  be  possible  to  assign  some  probabilities  to  one’s 
demand  forecasts  and  attempt  to  minimize  an  expected  present  value  function, 
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It  is  likely  to  prove  very  cumbersome  to  use.  In  practice,  planners  are 
more  likely  to  use  a  deterministic  model  and  treat  demand  forecasts  in  a 
parametric  fashion. 

There  usually  is  enough  time  between  expansions  to  permit  forecasts  to 
be  up-dated  and  new  determinations  made  of  optimal  policies  regardless  of 
which  model  is  used.  Whatever  the  length  at  the  planning  horizon,  if  demand 
is  not  a  linear  function,  a  stationary  policy  may  be  a  poor  approximation  to 
reality.  For  this  last  reason,  if  for  no  other,  we  feel  our  approach  yields 
results  which  may  be  of  practical  importance. 


CHAPTER  VI 


SOME  COMPUTATIONAL  RESULTS 


Introdurt Ion  , 

The  discrete  dynamic  programming  algorithm  for  Type  I  problems  was 
coded  in  Fortran  IV  for  the  Univac  1100  computer.  Computations  were  made 
for  seven  demand  functions. 

The  first  set  of  demand  functions,  shown  in  Figure  8  ,  were  chosen  to 
facilitate  comparing  a  stationary  policy  with  a  multi-stage  finite  planning 
interval  policy. 

The  second  set  of  demand-  functions,  shown  in  Figure  16,  were  chosen  to 

A 

test  the  sensitivity  of  x^  ,  the  optimal  value  for  the  first  stage  capacity 
expansion,  to  changes  in  demand  forecast,  using  a  family  of  demand  functions 
(all  of  which  coincide  during  the  first  third  of  the  planning  interval). 

The  following  were  chosen  so  as  to  simplify  comparisons  with  the  asymptotic 

case : 


a«.5,k«8,p*»l,r*0.1  ,1*60. 

All  computations  employed  the  same  length  planning  interval,  namely  60  units. 
The  units  of  measure  for  the  planning  interval  are  arbitrary.  The  demand 
functions  for  Cases  A,  B,  C  were:  linear,  strictly  concave  and  strictly  convex, 
with  D(0)  ■  0  and  D(I)  *  60  in  all  three  cases  see  Figure  8. 


Figure  9  shows  the  present  value  of  total  costs  as  a  function  of 


the  planning  interval.  The  value  of  n*  changes  at  20.  36,  and  51  periods 


Run  A:  Linear Demand  Function 


Present  Value  of 
Total  Cost  ($) 


FIGURE  9 
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Figure  10  is  the  lower  portion  of  Figure  9  and  is  drawn  to  the  sane 
scale.  It  shows  a  comparison  of  a  no-expansion  policy  with  the  policy  for 
n*  1  and  n  ■  2  expansions.  For  a  planning  interval  less  than  8.7  periods 
(approximately),  the  no-expansion  policy  is  optimal. 
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Case  A:  Linear  Demand  Function 


Comparison  of  No  Expansion  with  an  Optimal  n=l, 2  Expansion  Policy 


Present  Value  of 
Total  Cost  ($) 


Planning  Interval 


FIGURE  10 


Case  B:  Concave  Demand  Function 


Present  Value  of 
Total  Cost  ($) 


Planning  Interval 


FIGURE  11 


C a gp  f.  :  Convex  Demand  Function 


Present  Value  of 
Total  Cost  ($) 


FIGURE  12 
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Discussion  of  Computational  Results  for  Cases  A,  B  and  C 

The  computer  code  permits  a  maximum  of  seven  expansions.  In  Cases  A, 

* 

B  and  C,  the  optimal  number  of  expansions,  n  ,  was  four.  We  note  that 

the  optimal  value  of  the  cost  function  decreases  with  the  number  of  expansions 

*  * 
up  to  n  and  is  constant  for  all  n  >  n 

The  point  in  time  at  which  a  two  stage  policy  first  becomes  optimal  is 

almost  the  same  for  all  three  trials:  20  periods  for  concave  and  linear 

demand  functions,  22  periods  for  the  convex  demand  function.  Similarly, 

after  36  periods  a  three  stage  policy  is  optimal  for  Case  A,  B  and  C. 

Case  A:  Linear  Demand  Function 

In  Manne  [9],  Chapter  10,  Appendix  B,  the  following  results  are  given 
for  a  linear  demand  function  D(y)  «  y  using  the  asymptotic  stationary 
policy  cost  function: 

x*  -  15.17 
t*  -  3.12 

Y  -  4.287  . 

The  graphs  of  Figure  13,  showing  the  optimal  size  of  expansions  as  a 
function  of  the  length  of  the  planning  interval  illustrates  three  points  in 
time  for  which  the  optimal  policy  is  one  of  equal  expansions  —  30.34,  45.51 
and  60.68  periods.  These  points  coincide  with  the  "regeneration"  points  of 
the  stationary  model. 

Also  note  that  a  rather  long  planning  interval  must  be  used  before  a 
two-expansion  policy  is  optimal  (21  units)  and  that  smaller  planning 
intervals  may  imply  nonequally-sized  expansions. 
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For  planning  intervals  greater  than  21  periods  (when  the  two-stage 
policy  becomes  optimal) ,  the  optimal  size  of  the  first  expansion  varies 
only  slightly  from  the  stationary  value  of  15.17.  Similarly,  after  38 
periods  the  second  expansion  varies  only  slightly  from  its  stationary  value 
of  15.17.  Presumeably  similar  remarks  would  apply  for  the  third  expansion 
if  the  planning  interval  were  extended  beyond  60.68  periods.  Moreover,  the 
optimal  number  of  stages  for  60  periods  is  four,  the  same  as  for  the 
stationary  model.  These  findings  tend  to  confirm  a  principal  hypothesis  of 
a  stationary  model,  namely,  that  the  optimal  policy  is  one  in  which  the 
expansions  are  of  equal  size. 

Optimal  Size  of  Expansion 

In  this  section  we  consider  the  optimal  size  of  each  expansion  as  a 
function  of  the  planning  interval.  The  computational  results  are  in 
graphical  form.  (See  Figures  13,  14  and  15). 


i 
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Case  A:  Linear  Demand 


Optimal  Size  of  Expansion  as  a  Function  of  the  Planning  Interval  Length 
* 


FIGURE  13 

Since  only  discrete  values  were  computed,  the  continuous  curves  were  drawn  for 

^ 

ease  of  reading.  The  maximum  value  of  x^  ,  ,  x^  is  20.  The  stationary 

* 

value  of  x  ,  using  the  stationary  model, is  15.17. 

In  the  first  graph  above,  the  optimal  policy  is  x*  ■  0  for  a  planning 


interval  less  than  approximately  8.7  units  of  time. 
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Case  B:  Concave  Demand 
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Case  C:  Convex  Demand 

Optimal  Size  of  Expansion  as  a  Function  of  the  Planning  Interval  Length 


* 


* 


FIGURE  15 


Note : 

The  no-expansion  policy  was  not  computed;  hence,  graphs  are  drawn  as 

* 

if  the  optimal  policy  is  x^  >  0  even  for  small  planning  intervals. 


t 
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One  result  is  common  for  all  three  demand  functions,  namely  that 

A  A 

,  X2  the  sizes  of  the  first  and  second  expansions  tend  to  fluctuate 
around  some  constant  value.  And,  this  value  is  reached  at  a  point  in  time 

A 

just  after  the  optimal  number  of  stages  increases  from  one  to  two  for 

A 

and  from  two  to  three  for  x^  .  This  suggests  tha*  ,  for  example,  the  value 

A 

of  Xj  may  be  invariant  with  respect  to  a  family  of  demand  functions  which 
are  identical  rear  the  origin.  This  portion  is  usually  the  most  reliable 
part  of  a  demand  forecast.  Accordingly,  we  tested  this  idea  using  Case  C 
to  define  three  additional  demand  functions  denoted  Cases:  Cl,  C2  and  C3. 
These  four  demand  functions  are  identical  for  I  e  [0,20)  only.  (See 
Figure  16.)  These  functions  were  deliberately  chosen  to  provide  distinctively 
different  patterns  both  as  to  rate  of  growth  and  final  values  for  demand  at 
the  end  of  the  planning  interval. 


FIGURE  16 


SE 


STABILITY  OF  x} 


Cases:  C,  Cl,  C2  and  C3 


* 


FIGURE  17 


LEGEND 

-  Case  C 

-  Case  Cl 

- Case  C2 

—  .  —  Case  C3 


101 


The  results,  in  graphical  form,  are  shown  in  Figure  17.  Only  x^  ,  the 
optimal  value  for  the  first  expansion,  is  shown.  (For  our  purposes  we  have 
assumed  that  the  optimal  policy  is  not  a  no-expansion  policy.)  As  before,  a 

two  stage  policy  is  optimal  at  about  the  same  time,  between  22  and  25  years. 

*  * 

(This  coincides  with  the  first  decrease  in  x^  .)  In  each  case,  x^  tends  to 
a  constant  value  and  further,  these  constant  values  are  close  to  one  another. 

The  importance  of  this  result  is  as  follows.  In  practice  one  is  most 
interested  in  ascertaining  the  optimal  size  of  the  first  expansion  because 
it  is  the  first  commitment.  Subsequently,  new  demand  forecasts  and  a  new 
policy  can  be  determined.  It  is  generally  true  that  the  reliability  of  a 
demand  forecast  decreases  with  the  length  of  the  planning  interval.  This 
suggests  treating  demand  forecasts  parametrically.  Our  results  lead  us  to 
conclude  that  if  short  term  forecasts  are  reliable,  good  estimates  can  be 
made  for  the  optimal  size  of  the  first  expansion. 
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